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SUMMARY 


In order to provide a common basis for the comparison of high- 
performance supercharged aircraft engines, formulas are proposed 
for correcting the measured output to standard atmospheric con- 
ditions. 

An investigation of power-correction formulas is needed be- 
cause the formulas generally used were developed for naturally 
aspirated engines and have been found inaccurate when applied 
to supercharged-engine test data. 

In the development of the proposed correction formulas, the 
theoretic effects of atmospheric changes on the performance of an 
engine’s induction system units are discussed and applied in a 
general basic correction formula. These effects pertaining to each 
induction system unit are checked by comparing the unit formula 
with test data obtained on a specific engine. 

Two engine-output correction formulas, covering both full- 
throttle and part-throttle operation of single-stage supercharged 
engines, are then evolved. The accuracy of these formulas is 
determined by comparison of calculated engine performance with 
flight test data taken under various atmospheric conditions. A 
similar check of the correction formulas generally used is also 
made. 

A comparison of the old and new full-throttle formulas indi- 
cates that former correction errors of as much as 10 per cent may 
now be reduced to approximately 11/2 per cent. Comparison of 
the part-throttle correction formulas shows that the accuracy of 
the rational formula is slightly inferior to that of the generally 
used formula. 

Correction formulas for two-stage supercharged engines of the 
gear-driven and turbo-driven variety are finally proposed for con- 
sideration, but no flight test proof of these latter formulas is pre- 
sented. 

The conclusions reached as a result of the correction formula 
investigations are: 

1. The formula derived for use in correcting the full-throttle 
output of single-stage supercharged engines is basically more 
sound and considerably more accurate than the generally used 
full-throttle power correction formula. 

2. The rational formula derived for use in correcting the part- 
throttle output of single-stage supercharged engines shows no 
superiority over the empirical part-throttle correction formula 
now in general use. 
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3. The correction formulas proposed for use in connection 
with engines equipped with two-stage gear-driven and turbo- 
driven superchargers should give reasonable accuracy, although 
some experience must be gained in the use of these formulas 
before they are finally adopted. 


List OF SYMBOLS 


Cc = supercharger density parameter = 
V2 2.53 
cocoon er 
60887), 
C; = constant representing the ratio of the 


pressure head of cooling air on pass- 
ing through the intercooler, to the 





airplane pressure head, g. C, = 
| 52(p. — pa) 

‘ 1/9 (6«/32.17) V.? 

D (cu.in.) = engine displacement (swept volume) 

e = intercooler effectiveness = e = 
(Ta, — T)/(Ta, — Ts) 

g (ft. per sec.?) = acceleration due to gravity 

1 = subscript to denote ‘“‘observed”’ condi- 
tions 

k = ratio of specific heat values of air = 
C,/C. = 1.3947 

K. = equivalent orifice constant of carbure- 
tor at full-throttle setting 

K;. = equivalent orifice constant represent- 


ing the cooling air restriction losses 
through the intercooler and its cool- 
ing air ducts 

Ke = equivalent orifice constant represent- 
ing the engine air restriction losses 
through the intercooler and its 
engine air ducts 

N (r.p.m.) = engine crankshaft rotative speed 

n (r.p.m.) turbosupercharger impeller speed 


po (in. Hg) = absolute pressure in the atmosphere 

pd» (in. Hg) = absolute pressure at the main-stage 
supercharger entrance 

po, (in. Hg) = absolute pressure at the auxiliary- 
stage supercharger entrance 

p- (in. Hg) = absolute pressure at the carburetor 

: entrance 
pa (in. Hg) = absolute pressure at the main stage 


supercharger discharge annulus 
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pa, (in. Hg) 
pe (in. Hg) 
pr (in. Hg) 
dx (in. Hg) 
py (in. Hg) 
Qu (cu.ft. per sec.) 
Q, (cu.ft. per sec.) 
R, (ft.lbs. per Ib. 
per °F.) 


R, (ft.lbs. per Ib. 
per °F.) 


Sx (sq.ft.) 

S, (sq.ft.) 
FCB CE 
Ty (°R.) 

To, (“R.) 

T. (°R.) 

Ta (°R.) 


Ta, (°R.) 


Tx (°R.) 

IoC) 

V (ft. per sec.) 
V, (ft. per sec.) 
V, (ft. per sec.) 
V.. (ft. per sec.) 
V, (ft. per sec.) 
V,, (ft. per sec.) 


Wz (Ibs. per hour) 


W, (Ibs. per hour) 
W; (Ibs. per sec.) 


W; (lbs. per hour) 
W,, (Ibs. per sec.) 


Ww (Ibs. per hour) 
Y; 

Y, 

b.hp. (hp.) 

f.hp. (hp.) 

i.hp. (hp.) 


is.a.c. (Ibs. 
i.hp.-hour) 


per 
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absolute pressure at the auxiliary 
stage supercharger discharge annu- 
lus 

absolute exhaust back pressure 

absolute pressure at the cooling air 
exit face of the intercooler 

absolute pressure at the cooling air 
entrance face of the intercooler 

absolute pressure at the exit of the ram 
scoop 

volumetric flow rate of engine air 
through the intercooler 

volumetric flow rate of cooling air 
through the intercooler 


gas constant for air (53.35) 
gas constant for intake charge mixture 


subscript to denote ‘‘standard”’ condi- 
tions 

intercooler cooling air face area 

intercooler engine air face area 

absolute temperature of the atmos- 
phere 

absolute temperature at the main- 
stage supercharger entrance 

absolute temperature at the auxiliary- 
stage supercharger entrance 

absolute temperature at the carbure- 
tor entrance 

absolute temperature at the main- 

stage supercharger discharge annulus 

absolute temperature at the auxiliary- 
stage supercharger discharge annu- 
lus 

absolute temperature of the cooling 
air to the intercooler 

absolute temperature at the ram-scoop 
exit 

main-stage impeller tip velocity 

auxiliary-stage impeller tip velocity 

true airplane velocity 

velocity of the inducted air at the car- 
buretor 

velocity of the inducted air at the ram- 
scoop exit 

velocity of the inducted air at the 

intercooler engine air face 

mass flow of dry air inducted by en- 
gine 

flow of fuel to the engine 

mass flow of cooling air through the 
intercooler 

total flow of fluid inducted by engine 

mass flow of engine air through the 
intercooler 

mass flow of atmospheric water vapor 
into the engine 

main stage supercharger pressure 
parameter = [(pu/p,)**85 — 1] 

auxiliary stage supercharger pressure 
parameter = [(pu,/P»,)°°*83 — 1] 


ram scoop pressure parameter = 
[(pr/Pa)®*** — 1] 

engine brake horsepower output 

engine motoring friction horsepower 

engine indicated horsepower (b.hp. + 
f.hp.) 


indicated specific air consumption 





OF THE AERONAUTICAL SCIENCES—AUGUST, 





1942 






per cu.ft.) = density of inducted charge at entrance 
of main-stage impeller 

density of inducted charge at dis- 
charge annulus of main-stage super- 


6, (Ibs. 


6a (Ibs. per cu.ft.) = 


charger 

6, (lbs. per cu.ft.) = density of cooling air entering the 
intercooler cooling air face 

6, (Ibs. per cu.ft.) = density of inducted charge entering 
the intercooler engine air face 

n = main-stage supercharger pressure coef- 
ficient 

m = auxiliary-stage supercharger pressure 
coefficient 

Nr = ram-scoop pressure coefficient 

Ne = engine cylinder volumetric efficiency 


based on the density at the main- 
stage supercharger discharge annu- 


lus 

0 = symbol representing the ratio: 
(54/5) actuat/ (5a/5>)adiatatic 

ry = main stage supercharger temperature 
coefficient 

n’ = mechanical efficiency of engine at sea 
level 

r = ratio of mechanical friction to total 
friction horsepower at sea level 

1 = auxiliary-stage temperature coefficient 

INTRODUCTION 


AN A RESULT OF FLIGHT-TEST and test-stand investi- 
gations made some years ago on naturally aspi- 
rated engines,’ ? the following formulas were derived 
for the correction of engine output to standard atmos- 
pheric conditions. At full throttle 


pe, Fe ( i < mt) 
b.hp., = bhp. | — ./=-{1 + ——]} - 
‘ sie ke Ts n 


Cs] 0 


Because the ratio (A — An’/n’) had relatively minor 
importance at full-throttle conditions, it soon was al- 
most universally dropped. Hence, the full-throttle 
formula now in general use is: 


b.hp., = b-hp.;[(pes/pe:) V Tei/ Tes] (2) 


At part throttle, constant manifold pressure is as- 
sumed, and the term p,,/p,; does not apply. With con- 
stant pressure in the induction system, the generally 
used part-throttle correction formula becomes: 


b.hp., = b.hp.yW Te,/T, (3) 


Ever since the engine torque indicator has come into 
use in determining power output in flight, numerous 
cases have arisen in which the application of the full- 
throttle output-correction formula, Eq. (2), has defi- 
nitely given inaccurate corrected-output values. 

Because of this situation, an extensive investigation 
was conducted to discover the causes of power-correc- 
tion-formula inaccuracies. In this investigation both 
theoretic and practical studies were made and com- 
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It is the purpose of this paper to present the 


pared. 
results of these investigations. 


GENERAL ANALYSIS AND DEVELOPMENT OF THEORETIC 
POWER-CORRECTION FORMULAS 


In correcting engine output to standard conditions, 
two different points of view may be taken: The engine 
manufacturer generally wants to know what corrected 
output his engine will deliver if standard atmospheric 
pressure and temperature values are substituted for the 
observed values that existed at the induction-system 
entrance when the output was measured. This correc- 
tion of observed output to standard conditions provides 
a common basis on which engines of different design or 
manufacture can be compared. 

The airplane manufacturer is not especially interested 
in comparing engines. He wishes to know how his air- 
plane will perform in standard atmosphere but can sel- 
dom test it when the free air temperature, pressure and 
humidity are at arbitrary N.A.C.A. standard values. 
Therefore he must calculate the effect on airplane per- 
formance of correcting the engine output to correspond 
to standard atmospheric values at the airplane's air 
scoop entrance. This permits him to evaluate his air- 
plane’s performance with respect to the performance of 
other airplanes when flying in identical atmospheres. 

Both points of view will be considered in the develop- 
ment of the following engine-output correction formulas. 


Basis of Engine-Output Correction Formulas 


The fundamental basis of Otto-cycle engine-output 
relations is that which expresses the output as a func- 
tion of the air consumption. When the engine is oper- 
ated at fixed values of fuel-air-mixture ratio and spark 
setting, the engine power is found to vary in a linear 
relation to the dry mass airflow. That is: 

(4) 

In this relation the indicated specific air consumption 
is a constant that depends on the mixture-ratio, spark- 
advance setting and the combustion efficiency of the 
fuel mixture. As spark-setting and fuel-heating values 
are normally constant, it is primarily a function of the 
fuel-air-mixture ratio of the intake charge, although this 
function may be modified slightly by large changes of 
engine speed. 

The mass flow rate of fluid inducted by the engine 
cylinders is the sum of the flow rates of dry air, fuel and 


i.hp. = W,/i.s.a.c. 


the water-vapor content of the inducted air. This may 
be expressed by the relation 
W. Wy e 
W,=W.a\|1+ ek 1 + — | (approx.) (5) 
Wa Wa 


At the entrance of the engine intake manifold, the 
total fluid induction rate may be expressed by the 
following equation: 


W, = (D/1,728)(60N/2)n,54 (6) 
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Egs. (6), (5) and (4) may be combined, and the sum 
of the brake and friction horsepower may be substituted 
for the indicated horsepower to give the equation: 


b.ho. = — 0.01736D Nada 
ecsieie: 14 WT, 4 We 
1.8.a.C. - 
W, W, 


If the subscripts 7 and s are now used to denote observed 
and standard conditions of engine operation, respec- 
tively, the basic relation of standard to observed brake 
output is: 


b.hp., = [b.hp., + f.hp.,] X 


Ww, W 
N; 1 0 zt A.C. 4 l ~ 1 . 
(i (1+ “Al +H) 


Wo 
zt 

This is the formula on which the subsequently de- 
scribed methods of correcting engine output to standard 
conditions are based. It is a general equation that:is 
satisfactory for use with any four-stroke, Otto-cycle 
engine. 

In the development of correction formulas the prob- 
lem is to evaluate the different factors of Eq. (8) as 
functions of the atmospheric pressure, temperature and 
humidity so that the effects of atmospheric changes on 
engine output may be determined. 


—f.hp. (7) 


(8) 


— f.hp., 


W 
Namidajis.a.c.,( 1 + —)(1 4+ 
W.). 


DEVELOPMENT OF CORRECTION FACTORS FOR AN ENGINE 
WITH A SINGLE-STAGE SUPERCHARGER AND SUPER- 
CHARGER-ENTRANCE CARBURETOR 


For a conventional carbureted engine equipped with 
a single-stage supercharger, the engine-air induction 
and exhaust-gas disposal system is similar to that shown 
in Fig. 1. In the airplane installation, a forward-facing 
intake air scoop is generally employed to provide ‘‘ram,”’ 
as a result of the forward velocity of the airplane. A 
test-stand installation is more likely to provide an air- 
measuring and pressure-and-temperature regulating 
system than a forward-facing air scoop. 
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Fic. 1. Typical single-stage engine induction system diagram. 
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The basic difference between correction formulas for 
test-stand and airplane use lies in the fact that test- 
stand power corrections are made on the basis of the 
inducted air conditions at the entrance flange of the car- 
buretor, whereas the airplane power corrections are 
made on the basis of the free-air conditions at the en- 
trance of the forward-facing scoop. 


Full-Throttle Conditions 


At full throttle the inducted air undergoes numerous 
temperature and pressure changes before it reaches the 
entrance of the intake manifold. These changes are due 
to velocity variations in the air stream which are 
brought about by the supercharger and ram scoop, as 
well as by losses due to the carburetor chokes and 
throttles, the duct walls and elbows, and by the vapori- 
zation of the fuel as it mixes with the air stream. 

Changes in state of the inducted air due to these dif- 
ferent factors will now be considered. 


Scoop Effects on the Inducted Air 


The forward-facing airplane scoop produces both a 
pressure and a temperature rise which are functions of 
the scoop efficiency and the airplane speed, according to 
the relations :* 


| ( k—1V,? — \ 
r= pajl r - - i 
br = P. + \n b RT, 











be = 
V2 — V,2\?8 
Pa (1 + 1 ae) (9) 
T.=T,= rt + (A) "ee ) o 
k 2gR,T, 


Vi — V? 
ar) (10) 


T, (1 + 
Because the scoop efficiency is practically impossible 
to determine, except under actual test conditions, not 
only the relation of p, to p, but also the ratio 7, to T, 
must be measured during the flight test. Reasonably 
accurate means of estimating the carburetor-entrance 
pressure and temperature for standard conditions of the 
atmosphere are then available. For these estimations 
the scoop efficiency is assumed to remain constant de- 
spite the small variations in atmospheric conditions and 
airplane speed accompanying the correction proc- 
ess. 
For the airplane installation of a single-stage engine 
the relations of standard to observed carburetor- 
entrance pressures and temperatures are as follows: 


Pe, brs a | pe at mye) T" 
Pee Pt Pl + (Y,, ey =" 
hm me ee ee) SY 


‘aj? — Ve! 
<4 


a 
Va? — Ve? 


Fu 








Ta Tn te ~ Tap( lai 


Ti: Sa 
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These ram-scoop relations are not used in test-stand 
power corrections as the standard carburetor-entrance 
pressure and temperature are assumed to be equal to 
the standard free-air pressure and temperature values 
as determined from N.A.C.A. Report No. 218.4 


Carburetor Effects on the Inducted Air 


As the inducted air passes from the entrance scoop or 
duct on through the carburetor, there is a pressure loss 
created by the carburetor chokes and throttles. At 
part-throttle setting this loss is of large magnitude and, 
furthermore, is unpredictable. Fortunately, calcula- 
tions of the pressure drop through the carburetor need 
to be made only at full-throttle setting. A relation of 
the form 


Pe — Po = Kp. V.?/T, (13) 


will evaluate this full-throttle pressure drop. The 
ratio of standard to observed impeller-entrance pressure 
may be expressed in terms relative to conditions at the 
carburetor entrance as follows: 


pbs sd (2)(3 mer i # vel) 

pri a pa/\1 — K.VeP?/Te 
The temperature drop through the carburetor, due to 
fuel vaporization, is found to be dependent mainly upon 
the initial temperature of the entering air and the fuel- 
air ratio, regardless of throttle setting. Experimental 
data collected from several sources indicate that the 
drop in carburetor-air temperature due to vaporization 


may be calculated with sufficient accuracy by either of 
the empirical relations: 





(14) 


T, = 21.4VT, (15) 


or 


Ty = 0.521(t, + 873) (16) 


Other data (Fig. 2) indicate that changes in fuel-air- 
mixture ratio affect these calculated temperatures to 
only a minor degree and that further temperature drop 
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Effect of fuel-air ratio on the primary-stage impeller- 
entrance temperature. 


Fic. 2. 
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calculations involving the fuel-air ratio are unwar- 
ranted. Hence, the ratio of observed to standard im- 
peller-entrance temperature becomes: 


To;/Te, = WTei/Tes = (tee + 873)/(tee + 873) (17) 


Supercharger Effects on the Inducted Air 


On passing through the supercharger the inducted 
mixture of fuel and air is subjected to an appreciable 
pressure, temperature and density rise. When only 
air is the fluid worked upon, the accepted supercharger 
relations that express the pressure and temperature 
ratios as functions of the impeller tip speed and im- 
peller-entrance temperature are: 


Pa/P» = [(nV2/6,088T,) + 1]** (18) 


and 
Ta/T, = [(nV?/¢ 6,0887,) + 1] (19) 


Should perfect adiabatic compression be assumed, 
the pressure and temperature coefficients are equal to 
unity, by definition. Therefore, the adiabatic super- 
charger outlet-to-inlet density ratio is found by division 
of Eq. (18) by Eq. (19) as follows: 


lg 
dy adiabatic Po T4 adiabatic 


( V2 ) 
<n fe | 20 
6,088T7, 7 (20) 


Let the measured density rise ratio, divided by the 
calculated adiabatic density rise ratio, be a parameter 








represented by @. Then 
s (62/6»)actuar = (6a ‘Sy actuar —— 
(64 / by) adstadatic E tt he i. 
6,088 T, 


New 
b:)/ \Ta 
— eee (81) 
——__ + ] 
| emer, ” | 
Now let C represent the parameter [(V/?/6,0887,) + 
1)**. Then the supercharger-discharge density can be 


expressed in terms of the supercharger-entrance pres- 
sure, temperature and density by the equation 


= (2)(22) 
-, rhe 


ba = 6,0C 





or 
(23) 


For air, the pressure and temperature coefficients and 
density parameter are known to vary with the impeller load 
coefficient and may also be found to vary with initial tem- 
perature, especially at high impeller tip-speed values. 
However, for relatively small atmospheric pressure and 
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temperature variations the supercharger pressure and 
temperature coefficients and density parameter can nor- 
mally be assumed to remain constant. The ratio of stand- 
ard to observed supercharger-discharge or blower-rim 
pressures, temperatures and densities then may be 


written: 
me tts (myn YH (TA NM 
pas a 1! ni poi ViJ \T0,/. (24) 


Ta; 


i 


Ta, (Ta — To)(V./V;)? + Tr, 
5d,/bdg = (65,/50;)(C,/C;) 


(25) 


(26) 


By substitutions and combinations of Eqs. (14), 
(24), (25) and (26), the net theoretic effect of small 
changes in atmospheric conditions on the main-stage 
supercharger-discharge density, pressure and tempera- 
ture, may now be evaluated for the full-throttle condi- 
tions of the single-stage engine. 

For the test-stand corrections: 





Be (Ba) Fa\Y( GL = HVAT 
ba; a G i - seed ai 
me fs (E(B) 
pa; pai Ts J Ci pei 
Ve \2/ Te, \? 
ol 2Y(2)"T" es 
| + ( (4 Te; 
Ta; apes Ta; > 4 
Way BH gy vy, (29) 
(Ta, — 21.4Te*)\ 7") + 21.4Ts,"* 


The pressure and temperature at the forward-facing 
intake air-scoop exit are evaluated in separate equations 
so that Eqs. (27) to (29) inclusive, as well as the air- 
scoop relations, may be used for other more complicated 
conditions. For the airplane installation, the ratios of 
observed and standard values of the carburetor entrance 
pressures and temperatures required by Eqs. (27) to 
(29) are given by Eqs. (11) and (12). 


Intake Manifold Effect on the State of the Inducted Air 


In the intake manifold, the inducted charge may be 
subjected to still further changes in pressure and tem- 
perature. A well-designed manifold will be the cause 
of only minor pressure losses which, ordinarily, can be 
assumed to be negligible. However, because of adverse 
heat transfer losses, fluid friction losses, and residual- 
gas blow back, an appreciable temperature rise may be 
found at the intake port. Fortunately, this tempera- 
ture rise is generally found to be reproducible under 
any one set of operating conditions, although the magni- 
tude of the temperature increase may vary considerably 
with engine speed and, to some degree, with cylinder 
cooling-air flow and cooling-air temperature changes. 

Port-temperature variations with cooling conditions 
are almost impossible to predict with any accuracy. 
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Experience has shown, however, that this temperature 
function can be avoided by basing the cylinder pumping 
efficiency upon charge conditions at the main-stage 
supercharger outlet and by assuming that the same rela- 
tive increase of charge temperature in the manifold will 
occur under both observed and standard operating 
conditions. This procedure is believed to be adequate 
if the engine is operated at about the desired normal 
cylinder temperatures during the test period. Failure 
to do this may result in port-temperature variations of 
as much as 3 per cent. The effect is sometimes shown 
up by differences in corrected output which occur when 
flying both on a single engine and on two engines of a 
bimotored airplane. 


Volumetric Efficiency Variations with Atmospheric 

Changes 

Zeyns’ experiments’ have shown the engine cylinder 
volumetric efficiency to be a function of the engine 
speed, the valve timing, the intake-charge temperature 
and the ratio of the exhaust pressure to intake-manifold 
pressure. Of these different factors, the intake-charge 
temperature effect is apparently of least importance 
and can generally be neglected. The relation of super- 
charger-discharge volumetric efficiency to the above 
factors can be determined only by experiment on the 
type of engine under consideration. A typical set of 
volumetric efficiency curves for a conventional aircraft 
engine with normal valve timing is shown by Fig. 3. 
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Typical volumetric efficiency curves for a supercharged 
engine. 


Fic. 3. 


When operating under part-throttle conditions at 
constant pressure altitude, the ratio of exhaust back 
pressure can be expected to vary over a wide range 
of values, as the throttle opening and, hence, the mani- 
fold pressure is changed. This condition is typical of 
test-stand or sea-level dynamometer operation. Also, 
at constant manifold pressure and varying altitude 
conditions—i.e., varying exhaust back pressure— the 
manifold-to-back-pressure ratio is subject to rather 
wide variations. However, for a single-stage super- 
charged engine under full-throttle operation, the ratio 
of manifold pressure to exhaust back pressure is sub- 
stantially equal to the overall compression ratio of the 
supercharging system. This is, in turn, a function of 
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the impeller tip speed and the initial temperature at 
the supercharger entrance. As the ratio of exhaust 
pressure to supercharger-discharge pressure remains 
substantially constant for the full-throttle condition, 
the full-throttle volumetric efficiency of a single-stage 
engine, operating at fixed engine speed, can be assumed 
to remain substantially constant regardless of changes 
in atmospheric pressure or temperature. 

If it is desired to take the volumetric efficiency factor 
into account and if it is assumed that the observed and 
standard values of both the main stage supercharger- 
discharge pressure and the exhaust back pressure have 
been either measured or calculated, the corresponding 
observed and standard values of the cylinder volumetric 
efficiency may be found by reference to a curve similar 
to that in Fig. 3. The curve must be supplied by the 
engine manufacturer for the engine under consideration. 
On this curve the observed and standard values of volu- 
metric efficiency corresponding to the observed and 
standard values of the parameter (p,/p,) are spotted for 
the engine speed at which the power correction is being 
made. 

Recent experiments have indicated that the volu- 
metric efficiency may be varied as much as 2 to 4 per 
cent by wide changes in cylinder operating tempera- 
tures. Those changes are believed to be due largely to 
variations in tappet clearance as the cylinders expand 
orcontract. The effects of tappet clearance changes on 
output are unpredictable; therefore it is again suggested 
that attempts be made to test at the expected cylinder- 
operating temperatures, if possible, to eliminate errors 
due to cylinder temperature effects. 


Mixture Strength and Specific Air Consumption Variations 
with Atmospheric Changes 


In an aircraft engine equipped with a fully automatic 
carburetor, the indicated specific air consumption and 
the mixture strength are normally independent of at- 
mospheric conditions at the entrance of the induction 
system. The fuel-air-mixture ratio depends only upon 
the carburetor setting and the accuracy of the automatic 
altitude and economizer controls, and the specific air 
consumption is dependent primarily upon the mixture 
ratio. 

Unless the engine under test has been operated at 
other than the desired mixture-strength setting during 
the test period, there should ordinarily be no correction 
for fuel-air ratio or specific air consumption. However, 
under certain conditions of operation with automatic 
carburetors, this may not be the case. Fig. 6 illustrates 
a typical carburetor setting curve for an aircraft engine 
equipped with an automatically compensated carbu- 
retor. Should the engine be operating in either the “‘idle- 
transfer’ or ‘‘economizer’’ ranges under the conditions 
of test, then the correction to standard atmospheric 
conditions will involve not only a change in output, 
because of the change in inducted airflow, but also some 
change in output due to the change in mixture strength 
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which accompanies the change in airflow. It is appar- 
ent that not only the carburetor setting curve for the 
engine under consideration but also an estimate of the 
air consumption at test and standard conditions is re- 
quired for correction of output to standard atmosphere 
under these conditions. 

Inasmuch as specific air consumption and fuel-air 
ratio are closely interrelated, it is convenient to con- 
sider these two factors together rather than as separate 
entities. 

If all other factors of Eq. (6) are considered to remain 
constant, the variation of indicated horsepower with 
both mixture strength and specific air consumption be- 
comes 


if 1.8.a.C.; | (1 + (W,;/Wa)): 


~ is.a.c., [1 + (W,/Wa)s 


ihp.s 








ihp.; (30) 
Because the indicated specific air consumption must 
be determined experimentally anyway, it is convenient 
to determine by tests the empirical ratio of the i.hp. at 
any given mixture strength to the i.hp: at best-power 
mixture strength and to plot this ratio at the test mix- 
ture-strength value. This procedure implies that best- 
power mixture strength is the ‘‘standard”’ mixture. An 
experimentally determined curve of this sort is pre- 
sented to permit corrections to best-power mixture 
strength (see Fig. 5). 


Atmospheric Moisture-Content Effects 


Within the induction system the moisture content of 
the inducted air, in terms of the mass weight of water 
relative to the mass weight of air, can be expected to 
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remain constant, unless, of course, induction-system 
icing occurs. However, the density of the inducted 
charge may be appreciably affected by water-vapor 
content as well as by the supercharger, carburetor, etc. 

Under standard atmospheric conditions, based on 
N.A.C.A. standard air, the water-vapor content is 
assumed to be zero and the term [1 + (W,,/Wa)]s is 
equal to unity. The gas constant, which relates the 
inducted charge density to the pressure, temperature, 
water-vapor content and the fuel content of the charge 
may be determined by use of Fig.4. The standard gas- 
content value always lies on the extreme left line of this 
figure. 

The observed water-vapor content may be measured 
by means of such instrumentation as a hair hygrometer, 
wet-and-dry-bulb thermometer readings or a dew-point 
indicator. The ratio (W,/W,)i may then be found 


from any convenient psychometric chart with proper 
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Fic. 7. Variation of friction horsepower with engine speed, 
throttle setting and exhaust pressure. 


corrections for altitude effect. Assuming all other fac- 
tors constant, the correction for water vapor content 
may be expressed as follows: 


i.hp., Ww\ [Ru ‘ 
Ps (1 + *)( >) (31) 
i.hp.; Wa/ Ri: 

Engine-Speed Variations with Atmospheric Changes 


Practically all modern aircraft that are equipped with 
supercharged engines are also equipped with constant- 
speed propellers. The constant-speed propeller regu- 
lates the engine crank speed to the desired value regard- 
less of changes in atmospheric pressure temperature or 
density; therefore the ratio of standard to observed 
engine speed (V,/N;,) can ordinarily be assumed to be 
equal to unity. In some cases, however, it may be 
found that, because of oversight or a faulty tachometer, 
the engine was not operated at the desired crank speed 
and the effect of an engine-speed correction is desired. 
For this case the indicated engine-speed value is known 
and the standard engine-speed value is assumed. In 
making an engine-speed correction, it must be remem- 
bered that corresponding changes in supercharger tip 
speed must be included in the correction formulas. 
For gear-driven single-stage superchargers 


Vz, = ViAN;/N)) (32) 


and for the gear-driven auxiliary-stage supercharger, to 
be discussed later, 


Vig = Vi(Ns/ Ns) (33) 
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Stage Engine 


The friction horsepower of a single-stage supercharged 
aircraft engine, as determined by motoring tests on a 
dynamometer, is found to be dependent upon the engine 
crank speed, the oil and engine temperatures, the super- 
charger impeller-gear ratio, the airflow inducted by the 
engine and the exhaust back pressure. Fig. 7 is a typi- 
cal illustration of the effects of throttle setting and 
exhaust back pressure on friction horsepower. It is 
notable that the major effect of atmospheric changes on 
friction horsepower, all other factors being constant, is 
shown by the variation of motoring friction with exhaust 
back pressure. The variation of motoring friction with 
inducted airflow is found to be small with respect to the 
total friction and can be neglected for the small changes 
in inducted airflow which are represented by corrections 
from observed to standard atmospheric conditions. 
Minor oil or engine temperature changes are of negligible 
effect unless scavenging is also affected thereby. 

For all practical purposes, the relation of friction 
horsepower at standard conditions to that under ob- 
served conditions may be assumed to be a function of 
the exhaust back pressure only. The friction horse- 
power variation must be determined experimentally on 
a motoring dynamometer until more is known concern- 
ing the laws governing the effect of back pressure upon 
motoring friction. 


Single-Stage Engine Full-Throttle Power Correction 

Formulas 

By taking into account the effects of atmospheric 
changes on engine speed, friction horsepower, and the 
performance of various induction-system units, the 
theoretic full-throttle power-correction formula for a 
single-stage supercharged engine with constant-speed 
propeller may now be written as follows: 


N, Nos |[ Pes I’ i , 
b.hp., = [b.hp.; + thp.dl x [on [| 7 | ‘ 
reat 
HE BA 2d eed ( 3) x 


1 — (K.V-2/T.) || is.ac., |] /. . WA 


W; 
(2 7 a) 
W. Ru 
l — — —f.hp., (34 
i( . HAE i oT 


This relation may be simplified considerably by as- 
suming: (1) that the engine was tested at the desired 
engine speed—N,/N; = unity and V,/V; = unity; (2) 
that the exhaust-manifold back-pressure change is of 
small magnitude—»,,/n,, = unity and f-hp., = f-hp.i; 
(3) that the full-throttle loss through the carburetor is 
so small as to be of negligible effect, [1 — (K.V.,?/T.)/- 
[1 — (K.V.?2/T.,)] = unity; (4) that the mixture ratio 
is of the desired value and will not change appreciably 
because of economizer or idle-transfer action. 





C; 


Atmospheric Effects on Friction Horsepower of a Single- 























CORRECTION OF 
The above assumptions result in a simplified single- 
stage full-throttle formula, as follows: 


b.hp., = [b.-hp.; + f-hp.,] X 
We 
zm), 


bes (Tei\"* Cs 
Bs =) C; (1 bg 


For altitudes above 5,000 ft. the moisture content of 
the’air is generally low, so the terms 1 + (W,,/Wa) and 
R,/Ri may also be omitted. These terms are gen- 
erally much more important for dynamometer correc- 
tions than for flight-test corrections. 

When the engine is installed in an airplane, the values 
of p../P, and T,,/T., are determined by Eqs. (11) and 
(12). For convenience, the chart of Fig. 8 is presented 
for the determination of the parameter C. 


Ru 


Zi —f.hp.; (35) 
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Fic. 8. Chart for rapid determination of the parameter ‘“‘C.”’ 


Part-Throttle Single-Stage Correction Formulas 


At part throttle, for which condition the intake- 
manifold pressure is assumed to remain constant, a 
theoretic power-correction formula can be written as 
follows: 


N, os Ta W[i-s.a.c. 
nip, = bain + top [S| A Be] 
i vi 8 ‘silt laae 
W; 
1 + 
( =H We Ri 
peo. —|—fhp., (36 
L(+ wi). Le] - e+ 


st 
(1 + 7) 
‘d/s 
For substitution in this formula, the ratio Ty,/T,, has 
been evaluated in terms of 7., and T,, by Eq. (29). 
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ENGINE OUTPUT 


If the same simplifying assumptions are made for the 
part-throttle condition that were made for the full- 
throttle condition, the theoretic part-throttle correction 
formula may be reduced to the simplified form: 


b.hp., = 


Te | 
eee 
aise th Miia Bs - eT. = T. ‘J ‘o 


1+ =) Bee es 37) 
L(+ ea) Lee) = e+ 


It will be shown later that part-throttle Eq. (37) is 
correct in form but that the temperature term is incor- 
rect because of other factors not taken into account in 
the derivation of the equation. 


ay a ony Checks of the Single-Stage Power Correction 
‘ormulas 

In order to check the rationally developed single- 
stage power correction formulas, two special test setups 
were made. An engine was first thoroughly calibrated 
on an electric cradle dynamometer stand after which it 
was installed in a Vultee AB-2 airplane. The power 
plant used for the investigation was a Pratt & Whitney 
R-1830 engine equipped with a single-stage two-speed 
supercharger. 

Besides the usual calibrated engine instruments, 
several items of special equipment were installed to per- 
mit the measurement and control of the following: 


Item Instrumentation Control 
Engine speed Sensitive tachometer (d) Electric dyna- 
mometer 
(a) Constant speed 
propeller 
Engine output (d) Cradle dyna- 
mometer 
(a) Torquemeter and Throttle, etc. 
gauge 
Induction system (d) Mercury man- (d) Air chute throttle 
pressure ometers 


(a) Sensitive manifold (a) Altitude change 
pressure gauge 
Induction system (a) Thermocouples 
temperatures & and 
(d) Potentiometer 
(d) Mercury man- (a) Throttle in ex- 
ometer &  haust manifold 
(a) Sensitive manifold (d) tail pipe 
pressure gauge 
(a) Calibrated orifice 


(a) Carburetor  pre- 
heat system 


Exhaust system 


pressure 


Induction system (a) Throttle and den- 


airflow & and & sity change (al- 
(d) pressure gauge (d)  titude and pre- 
heat) 
Cylinder tem- (a) Thermocouples (d) Cooling blower 
peratures & and throttle 


(d) Potentiometer (a) Cowl flaps 
Nore: (a) represents airplane installation and (d) represents 
dynamometer installation. 
The airplane test setup is illustrated by Fig. 9. 
On the dynamometer, numerous so-called “‘altitude’’ 
calibrations were made for varying conditions of exhaust 


back pressure, mixture ratio, cylinder cooling, carbu- 
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retor-entrance pressure, impeller-gear ratio and engine 
speed. Because of unsatisfactory carburetor-air tem- 
perature control, the dynamometer tests were all made 
at a constant carburetor-air temperature of 90°F. 

The flight tests consisted of full-throttle runs at 1,800 


and 2,300 engine r.p.m. at altitudes of 5,000, 10,000, 
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Fic. 10. Typical performance and output-correction curves 


obtained during temperature-test investigations. 
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Flight test setup for power correction-factor investigations. 


14,000 and 18,000 ft. At each altitude the supercharger 
was operated in both low and high impeller-gear ratios. 
At 5,000 and 10,000-ft. altitudes, the variation of engine 
output with carburetor-air temperature was determined 
by varying the carburetor preheat valve from full-cold 
to full-hot positions, the range of temperatures being 
of the order of 60° to 80°F. at any one altitude. For 
these runs, the exhaust-manifold throttle valve, as well 
as the carburetor throttle, was locked wide open. 

Another series of runs at 5,000 and 10,000 ft. was 
made with the preheat valve in the full-cold position 
and with the exhaust-manifold throttle varied to pro- 
duce back-pressure changes of as much as 15 in. Hg. 

For all tests, the engine cylinder temperatures were 
kept at about the normal value for continuous opera- 
tion (450°F.) by control of cow! flaps. The carburetor 
mixture was controlled manually to obtain ‘‘best power’’ 
setting at each operating point. 

At each different test condition, the induction system 
pressures, and temperatures, the cylinder temperatures 
and baffle-pressure head, the engine speed and output, 
the exhaust-system pressure, the inducted airflow, fuel 
flow and water-vapor content were measured and re- 
corded. 

The test data taken for varying carburetor-air tem- 
peratures were plotted as shown by the typical curves 
of Fig. 10. Corrections to standard pressure and tem- 
perature at the carburetor entrance were made to the 
observed engine output values, using the precise Eq. 
(34), the simplified Eq. (35) and the generally-used Eq. 
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(2). The manifold pressure was also corrected to 
standard conditions by use of a simplified form of Eq. 
(28)—namely, 


Pas = Pea {1 + [Vi (Ta/Tes)”] f°” 


The results of these corrections are also shown on Fig. 
10. The superiority of the proposed full-throttle Eqs. 
(34) and (35) to the generally used Eq. (2) is to be noted. 
Typical results of the back-pressure runs are shown 
by Fig. 11. Each point on this curve was corrected to 
standard atmospheric back-pressure conditions by 
correcting, first, for volumetric efficiency variations 
with p,/pa and, second, for friction horsepower varia- 
tions with p,. The experimental volumetric efficiency 
and friction-horsepower data used for these corrections 
were determined from the extensive test-stand runs on 
the engine (Figs. 3 and 7). The excellence of the back- 
pressure correction method is also to be noted. 


(38) 
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Fic. 11. Typical performance and output-correction curves 


obtained during back-pressure test investigations. 


In order to check the assumptions that the super- 
charger pressure and temperature coefficients and the 
density parameter remain sufficiently constant during 
power correction periods to justify the correction 
formulas, the experimentally determined values were 
plotted against the carburetor-entrance temperature, 
as shown by Fig. 12. It is notable that for any fixed 
engine speed the pressure coefficient and density parame- 
ter remained substantially constant regardless of 
carburetor air temperature. However, the tempera- 
ture coefficient was found to vary over a wide range of 
values, some of which were unbelievable. 

To determine the cause of this extreme variation of 
temperature coefficient with carburetor-entrance tem- 
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perature, all induction-system temperatures were 
plotted against carburetor-entrance temperature as 
shown by Fig. 13. Inspection of these plots led to the 
conclusion that for low carburetor-entrance tempera- 
ture values the fuel continues to vaporize as it passes 
through the impeller and diffuser. * This gives an effec- 
tive temperature drop, below the carburetor air tem- 
perature, very much as indicated by the lower dashed 
line in Fig. 13. 

For best-power mixture, this effective drop amounts 
to about 70°F. at an initial carburetor air temperature 
of 130°F., and 45°F. at a carburetor air temperature 
value of 0°F. For the engine model tested, a formula. 


of the form 


T, = 0.812(T, + 50) (39) 


expresses the relation of the effective supercharger- 
entrance temperature with respect to the carburetor- 
entrance temperature. 

The relation of observed to standard effective im- 
peller-entrance temperature for the above relation be- 
comes: 


Ti Ts ” (Te, 50) (7. + 50) (40) 
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and a revised formula for the supercharger-discharge 
temperature at constant impeller tip speed is as follows: 


Tai/Tas = Ta/(Ta + 0.812(T., — Tei)] (41) 


This relation, instead of Eq. (29), when combined with 
Eq. (36) and simplified as before, should be expected 
to give a more accurate part-throttle correction formula 
to replace Eq. (37) as follows: 


b.hp., = [b.hp.; + f-hp.,] X 


To 1 We) [Re] _ pp 
s + 0.812(T2, — Ta) | Wa te | want 


(42) 





PART-THROTTLE CORRECTION FORMULA CHECKS 


In order to check the part-throttle, constant mani- 
fold pressure, formulas by use of the full-throttle test 
data of Fig. 10, the engine output at the different carbu- 
retor air temperatures of this figure were first corrected 
to constant manifold pressure by the relation: 


b.Hp.(corr.) = (b.hp. ; + f.hp.;)(Pas/Pai) = f.hp., (43) 
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For this relation the value of pz, was assumed to be 39.8 


in. Hg. A plot of these pressure-corrected output val- 
ues is shown on Fig. 14. The part-throttle correction 
formulas were next applied to the above-plotted output 
values. These results are also shown on Fig. 14 and 
the relative accuracy of the various formulas is demon- 
strated. 

It is apparent that neither of the rational part-throttle 
formulas gave as good accuracy as the empirical Eq. (:). 
The three part-throttle formulas were also applied to 
test data taken on other engine models. In practically 
all cases, empirical Eq. (3) was as good as or better than 
either Eq. (37) or Eq. (42). 


Check of Assumptions in Derivation of Formulas 


Several other plots of the test data were made to 
check some of the basic assumptions made in the deriva- 
tion of the correction formulas. In order to check Eq. 
(4), all the indicated horsepower values obtained at 
best-power mixture ratio were plotted against the dry 
mass airflow as determined by the calibrated orifices. 
The results are given in Fig. 15. The direct linear rela- 
tion between the indicated horsepower and the dry mass 
airflow justifies the use of Eq. (4). 





1400 
































= im 
af. Ba vd 
| | | | 
3000 4000 5000 6000 7000 8000 9000 

DRY AIR FLOW, Wg LB/HR. 











Linear variation of indicated horsepower with dry 
mass airflow. 


Fic. 15. 


The relation of indicated specific air consumption to 
fuel-air ratio, engine speed, output, etc., was found by 
plotting several test i.s.a.c. values against mixture ratio 
for the various engine-speed and engine-output values 
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tested. This plot, Fig. 16, indicates that the mixture 
ratio is the major variable affecting the indicated specific 
air consumption. Variations of i.s.a.c. due to engine 
speed and engine output were practically within the ex- 
perimental error of testing. 
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Variation of indicated specific air consumption with 
fuel-air mixture ratio. 
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The curves of Figs. 2, 3, 5 and 7 were also the results 
of plotting the test data obtained during the dyna- 
mometer and flight tests. 

Because of the generally satisfactory correlation be- 
tween the theoretic formulas and the experimental re- 
sults obtained on the single-stage engine, theoretic 
studies of correction formulas for two-stage engines 
with gear-driven and turbo-driven auxiliary stages were 
continued. Unfortunately, no carefully controlled 
flight-test checks of these formulas were possible. 


CORRECTION FORMULAS FOR THE TWO-STAGE SUPER- 
CHARGED ENGINE 


The induction system of a typical airplane installa- 
tion of an engine, equipped with a two-stage gear- 
driven supercharger, is illustrated in Fig. 17. Between 
the carburetor and the intake air scoop of the single- 
stage engine another supercharger and an intake-charge 
intercooler have been added. At full-throttle condi- 
tions, the ram-scoop-exit pressure and temperature are 
substantially equal to the pressure and temperature at 
the entrance of the auxiliary-stage impeller. That is, 
bn = p,and T;, = T,;. Eqs. (9) and (10) again may be 
used to calculate the values of 7, and #,. 

The pressure and temperature at the auxiliary-stage 
outlet can be determined by relations similar to those of 
Eqs. (18) and (19): 





Typical induction system of a two-stage engine with 
direct gear-driven auxiliary stage. 


(44) 


; . V,2 3.53 
Pa, fe (OEE 
Po, pr 6,0887, 


Ta _ Ta 


T,, ~ T, (45) 


mV? 
‘ws € 6,088T, * 1) 

For the auxiliary-stage supercharger the use of a 
density parameter @ is unwarranted, since the pressures 
and temperatures must be considered as separate enti- 
ties in correcting to standard conditions. 

On passing through the intercooler, the air dis- 
charged from the auxiliary-stage supercharger under- 
goes both a pressure and temperature change. The 
pressure change is largely dependent upon the air den- 
sity and the mean velocity at the induction-air face of 
the intercooler, whereas the temperature drop is de- 
pendent upon the temperature and mass flow of cooling 
air passing through the intercooler in the cross-flow 
direction. 

The pressure drop of the inducted air passing through 
the intercooler may be determined by a relation similar 
to that of Eq. (13) 


Pa ia P- = Kupa Vu?/Ta, (46) 


In estimating the variation of the engine air tempera- 
ture drop through the intercooler due to atmospheric 
changes, it is necessary to make some assumptions con- 
cerning the change in cooling air conditions in the inter- 
cooler as well as the change in engine air conditions. 

At a fixed ratio of mass cooling airflow to mass engine 
airflow, the intercooler effectiveness, e, is given by the 
relation: 

(47 


e= (Ta, — T.)/(Ta — Tr) 


This intercooler effectiveness can be expected to re- 
main substantially constant over a reasonably wide 
range of pressure and temperature conditions if the 
ratio of inducted engine airflow to cooling airflow does 
not change appreciably. 

For relatively small changes of atmospheric pressure 
and temperature, tiuus latter condition will be found to 
be approximately true if the overall conductivity of the 
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intercooler and its cooling air duct is not changed by the 
movement of the cooling-duct exit-gill-conirol flap and if 
the airplane speed does not change appreciably due to 
the expected change in engine output. 

The mass flow in the cooling air direction is given by 
the formula: 


W, = O15, saad C.K Sy V6; (48) 


and the mass flow in the engine air direction is given by 
the formula: 


: 141.6(pa, — po) 
W. = Oud, = ae eo -KySyb, (49) 
From Eq. (46) 
Pa — Pe = 0.7546,K, V,,? (50) 
Therefore, 
W, = 58.5K,'"S, Vib, (51) 
and 
of = Soren: = Const. X (52) 


At constant engine speed and full-throttle setting, the 
mean velocity of the air at the intercooler engine air 
face (V,) will be found to be relatively independent of 
atmospheric density changes in the induction system, 
and both 6, and 6, will vary approximately in proportion 
to the atmospheric air density 5,. If the airplane veloc- 
ity change is also assumed to be negligible, it is appar- 
ent that the ratio of cooling airflow to inducted engine 
airflow is approximately constant and the function e 
may also be assumed to be constant for correction for- 
mula derivations. 

If Eq. (47) is now solved for 7, the result is: 


T, = Ta,(1 — e) + eT, (53) 


Inspection of Eqs. (12) and (44) will show that for con- 
stant auxiliary-stage impeller tip speed and for the 
small changes in airplane speed and scoop-exit velocity 
ordinarily involved in correcting engine output to 
standard conditions the temperature rise in the ram 
scoop and in the auxiliary-stage supercharger remains 
constant regardless of initial temperature or pressure. 
Similar reasoning will show that the temperature rise 
of the intercooler cooling air, as it enters the cooling air- 
entrance duct, also remains substantially constant re- 
gardless of initial air temperature or pressure. 

These statements may be expressed algebraically as 
follows: 

Ta, = Const.; + Const. + 7, (54) 

and 


(55) 


The substitution of these values in Eq. (53) gives the 
following equation: 


T, = Const.; + 7, 
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T, = (Const.; + Const.,)(1 — e) + 
(1 — e)T, + e Const.; + eT, 


(56) 
or 
T. = T, + Const.; (57) 

Eq. (57) shows that for the purposes of correcting 
engine output to standard conditions the net tempera- 
ture rise of the atmospheric air in passing through the 
ram scoop, auxiliary-stage supercharger and intercooler 
can be assumed to be constant at fixed engine-speed 
condition. 

For constant engine speed the relation of observed 
to standard carburetor air temperature may now be 
written: 


To te = T/ (Teg 7 Ti + Tas) (58) 


Also, the standard temperatures in other parts of the 
auxiliary-stage induction system may be calculated by 
the relations 


Ta, = Tay, — (59) 


Ta + T., 
and 


T,, wei T,, = Ta; i Tas (60) 


These latter temperatures must be used in determining 
the standard pressures at various locations in the induc- 
tion system. The standard pressure at the auxiliary- 
stage discharge is given in terms of the ram-scoop 
pressure and temperature by a relation similar to Eq. 
(24): 


Pars = Proll + (Yi Vig/Vis)*(T/T 1.) ]°™ 


The standard pressure at the carburetor entrance may 
also be given in terms of the auxiliary-stage discharge 
pressure and temperature by the relation: 


re _ (bau = ba\(Tan\(Vua\*] (G9 
Pes = Péis E ( par; Nehe) | i: 


The simplifying assumptions that changes in airplane 
speed, scoop-exit velocity and induction-air velocity at 
the intercooler engine air face are negligible may now be 
made, and Eqs. (11), (12), (58), (59), (60), (61) and 
(62) may be combined to give the relation 


[i+ Ye (rt) [x 

’ Fi 73.53 
li! i (F. — Ty + ma ‘ 
= (Gee: 
pa; Ta; 


(61) 





Pes _ Pas 


Pei Pei 


= Tas + ome (69) 


This equation may be combined with Eqs. (58) and (55) 
to obtain a somewhat complicated full-throttle correc- 
tion formula for airplane-installed engines equipped 
with gear-driven two-stage superchargers: 
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(64) 





F= 


It is notable that it is necessary to know the test 
values of practically all the induction-system pressures 
and temperatures in order to correct the two-stage 
engine output to standard conditions. 

Eq. (64) covers the case in which the cooling air pass- 
ing through the two-stage engine intercooler is unregu- 
lated by intercooler cooling air control flaps. If either 
automatic or manual control of the cooling air to the 
intercooler is assumed, it is reasonable to suppose that 
such control will be for the purpose of obtaining con- 
stant carburetor air temperature. When carburetor 
air temperature control is assumed, atmospheric pres- 
sure and temperature changes will affect only the car- 
buretor-entrance pressure. In order to correct the two- 
stage full-throttle engine output to standard conditions 
under this type of operation, the terms [7.,/(T., — 
Ty, + To)|'* and C,/C; of Eq. (64) should be 
dropped. 

A full-throttle correction formula for a dynamometer 
installation of a two-stage engine has not been derived, 
since the two-stage engine is ordinarily supercharged so 
highly that full-throttle operation at or near sea level 
is impossible because of engine strength and detonation 
limitations. 

When a two-stage engine is tested on a dyna- 
mometer with the auxiliary stage in operation, the 
air-throttle and carburetor-entrance temperatures are 
either regulated to standard atmospheric pressure and 
temperature conditions or the carburetor-entrance pres- 
sure is regulated to the desired value and corrections are 
made only for carburetor air temperature and back- 
pressure variations from standard. These corrections 
simulate the so-called part-throttle corrections for the 
flight installation. 


PART-THROTTLE CORRECTION FACTORS FOR GEAR- 
DRIVEN TWwo-STAGE ENGINES IN FLIGHT 


The two-stage engine, when operating below its 
maximum critical altitude, may be in “neutral,’’ “low 
auxiliary” or “‘high auxiliary’’ condition. 

When the engine is operating in the neutral condition, 
it is essentially comparable to a single-stage engine, with 
the exception that for the neutral condition scoop 
losses, from atmospheric air pressure to carburetor air 
pressure, are likely to be relatively large. For flights 
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up to the main-stage critical altitude, with the auxiliary 
stage in neutral, the correction Eq. (42) applies. Above 
the main-stage critical altitude Eq. (35) applies until 
the auxiliary stage is shifted into the low-auxiliary 
condition. 

Immediately after this shift is made, 
matic carburetor-entrance pressure control is brought 
into operation, so that below the low-auxiliary 
critical altitude p,, = »,, and the part-throttle correc- 
tion formula becomes: 


b.kp.. = [bh |) ot np 
P- aed P- ee ere “ 


C, Ww\ (R: 
(E)C + W “\(R) | —f.hp.; (65) 


For this case the assumption that the temperature 
rise from the outside air to the carburetor entrance is 
constant is probably not quite so justified as it is for the 
full-throttle condition. For the part-throttle condition 
the ratio of mass airflows through the two respective 
entrance faces of the intercooler probably does not: re- 
main so constant as it does for the full (air) throttle 
condition. 

On reaching the full-throttle condition in low-auxili- 
ary stage the automatic carburetor-entrance pressure 
control becomes inoperative and Eq. (64) applies. 

When the auxiliary-stage impeller is next shifted into 
high-auxiliary position, the automatic carburetor- 
entrance pressure control again becomes operative and 
Eq. (65) applies up to the full air-throttle opening posi- 
tion, when Eq. (64) again applies. 

In using the above-mentioned correction formulas, it 
must be remembered that the observed friction horse- 
power, f.hp.;, and the auxiliary-stage impeller tip-speed 
change to a different value at each different auxiliary- 
supercharger operating condition. 


an auto- 





Correction Formulas for Turbosupercharged Engines 


Fig. 18 illustrates a typical induction-system dia- 
gram for a two-stage engine with a turbo-driven auxili- 

























Typical induction system of a two-stage engine with 
turbo-driven auxiliary stage. 


Fic. 18. 
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ary-stage supercharger. As far as correction formulas 
are concerned the basic difference between the turbo- 
driven and the gear-driven superchargers lies in the fact 
that the tip speed of the turbosupercharger impeller 
does not necessarily remain constant, even though the 
engine speed may be held constant. When the turbo 
is in operation below critical altitude, the automatic 
turbo-control is ordinarily adjusted to maintain a con- 
stant manifold pressure by varying the turbo-exhaust 
gas waste gate position. Critical altitude for the 
turboengine generally occurs when maximum allowable 
turbowheel tip speed is reached. It may be assumed 
that above critical altitude the turbowheel and turbo- 
supercharger impeller tip speed is regulated to a sub- 
Stantially constant limiting value and that Eq. (64) ap- 
plies in output corrections. 

In this case constant engine-speed operation is exactly 
similar to gear-driven auxiliary supercharger operation, 
with the exception that the turbosupercharger tip- 
speed value is fixed regardless of the choice of engine 
speed. 

A correction formula that applies to operation of a 
turbo-two-stage engine at altitudes below the turbo- 
critical altitude or critical turbo-speed is exceedingly 
complex and is therefore too cumbersome to use. 
This is because the operation and control of the turbo- 
stage is quite complicated. In attempting the deriva- 
tion of the correction formula for ‘‘below-critical-alti- 
tude’’ use, it was assumed that the turbocontrol would 
be such as to maintain constant manifold pressure at the 
main stage supercharger-discharge annulus. For this 
reason, the correction will be termed the ‘‘constant- 
manifold-pressure’’ correction. 

In order to simplify the correction formula as much 
as possible, it was also assumed that all corrections 
should be made at constant pressure altitude—i.e., 
Pas = Pa. The constant manifold-pressure output cor- 
rection formula, therefore, concerns the effects of 
changes of atmospheric air temperature on the turbo- 
engine output. 

A brief explanation of the turbo-operation for changes 
in free-air temperature is deemed advisable to show its 
complexity. 

Let it be assumed that the engine was tested 
in an atmospheric temperature somewhat above 
the standard value at the test pressure altitude. For 
this case the temperature at the entrance of both the 
primary and auxiliary stages would be in excess of the 
standard value, and the corresponding pressure rise 
ratio would be low in each supercharger unit. In order 
to obtain a definite manifold pressure, the turbowheel 
speed would be adjusted to a value somewhat higher 
than would be required if the atmospheric temperature 
were at its standard value. 

A complicated relation for calculating the ratio of the 
turbowheel speeds under observed and standard condi- 
tions is as follows: 
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(approx.) 
An unpredictable exhaust back-pressure change also 
accompanies the turbowheel speed change. 
The*approximate standard temperatures within the 


induction system at varying locations are given by 
Eqs. (29), (60), and the following equations: 











(66) 
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Ta, = (Ta; = T)( (67) 





(68) 


The combination of these relations gives an approxi- 
mate formula for calculating the standard value of 
main-stage supercharger-discharge temperature: 


Vi,\2 
Ta = Ta — 0.812T., + o.s12| ( (Ta ~t (3) “ 


lj 
Ta ~— Te 
I at Fa y am — aa 
+i 


Ta, — y 4 1/s 
kj Gey (approx.) (69) 


In the preceding equation the ratio V;,/V;,; is given by 
Eq. (66). In order to obtain an actual final form of an 
equation for determining the standard value of main- 
stage supercharger-discharge temperature, the solution 
of some extremely complicated simultaneous equations 
isrequired. The relation of standard to observed horse- 
power would then be given by the relation 


b.hp., = (b.hp.; + f.hp.,)(Tx/Ta,) — f-hp.s 


In view of the requirements of numerous assumptions 
and the complex nature of the derived basic equations, 
which involve many test measurements and unpredict- 
able factors, it is believed that the presentation of a final 
formula for the constant manifold-pressure correction 
of turboengine output is unwarranted. Hence, none 
is included in this paper. 


(70) 


CONCLUSIONS 


As a result of the analytic and experimental investiga- 
tions of the problem of correcting engine output to 
standard conditions, it is concluded: 

(1) That the commonly used full-throttle power 
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correction formula does not give accurate answers when 
applied to supercharged engine test data. 

(2) That the commonly used part-throttle correc- 
tion formula gives reasonably accurate answers when 
applied to single-stage engine test data. 

(3) A rational full-throttle output correction for- 
mula for use with single-stage supercharged engines has 
been evolved and is presented as Eq. (35) in this paper. 
Experimental checks of this formula indicate that its 
accuracy is considerabjy better than the accuracy of the 
commonly used full-throttle power correction for- 
mula. 

(4) A rational part-throttle output correction for- 
mula for use with single-stage supercharged engines has 
been evolved and is presented as Eq. (42) in this paper. 
Experimental checks of this formula indicate that its 
accuracy is inferior to that of the commonly used part- 
throttle correction formula. 

(5) Rational full-throttle and part-throttle correc- 
tion formulas are presented for use with two-stage 
engines of both the gear-driven and turbo-driven vari- 
ety. For the direct gear-driven two-stage engine, Eqs. 
(64) and (65) are expected to give reasonably accurate 
results for full-throttle and part-throttle corrections, 


Dear Sir: 

The letter from W. K. Bodger in the March JouRNAL, 
commenting on the paper, ‘‘Determination of Strain 
Distribution by the Photo-Grid Process,”’ is particularly 
gratifying to the authors, as it indicates an interest in 
the subject at an outstanding educational institution. 
It frequently occurs in industry that certain lines of 
research are limited to practical applications in that 
specific industry. When this point is reached it is of 
value to industry in general if the subject can be carried 
forward by a university or an institution having facili- 
ties for long-range research. The photo-grid process is 
being extensively used at Lockheed as a method of 
strain determination in a variety of forming operations, 
as described in ‘‘Elastic Theory Applied to Sheet Metal 
Forming,’ by F. R. Shanley.!. However, a number of 
long-range projects utilizing the photo-grid system to 
study local deformation adjacent to cutouts and notches 
in both impact and static test specimens would be of 
value. Utilizing high-speed light sources it should be 
possible to measure the strain distribution over impact 
specimens similar to the diagram obtained at normal 
speed as shown in Fig. 1 of this paper, using a specimen 
similar to Fig. 2, reference 2. 

Using the photo-grid process, finer grids may be 
placed over areas to be investigated in a fraction of the 
time required by a ruling engine. Furthermore, the 
deleterious effect of the scribed notch is not present. 
It would seem that the proper place for investigation 
into subjects similar to those suggested would be at a 
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although no carefully controlled experiments are pre- 
sented as proof of these formulas. 

6. Output corrections for constant manifold-pres- 
sure operation of a turbosupercharged engine are so 
involved as to preclude the possibility of any accuracy 
of correction. 

7. The full-throttle output-correction formula (Eq. 
64) derived for the two-stage engine with gear-driven 
superchargers applies also to turbosupercharged en- 
gines, at altitudes above the critical altitude, where 
constant turbospeed is assumed. 
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Enlarged approximately 38 diameters. 


Fic. 1. 
university, since the aircraft industry is at present 
preoccupied with the solution of immediate practical 


problems. 
Mr. Bodger has pointed out that true strain has been 


defined as 
L 
«= I dL/L = 1n(L/L») 
L 


0 


where Ly equals the initial gauge length and L equals 
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Fic. 2. Enlarged approximately 90 diameters. 


the final gauge length (reference 3). In the same 


reference it has been shown that 
e = In(L/Lo) = In(A/Ao) 


Therefore, if the initial and final areas of the cross 
section are known, the true strain at that section may 
be calculated from the relationships given above. Mr. 
Bodger has suggested that the “‘hypothetical ultimate”’ 
strain referred to in our paper (reference 2, p. 4) was 
indeed the true strain at fracture; and that whereas we 
stated that this strain could only be determined by 
using a gauge length of differential size, actually this 
strain could be measured by the substitution of the 
measured initial and final areas into the relationship 
In(A/Ao). Unfortunately, aircraft sheet mate- 
rials necessarily involve material thicknesses that pre- 
clude the possibility of making round tensile specimens 
from the sheet, so that flat coupons following A.S.T.M. 
practice must be used. Utilizing the flat specimen it is 
difficult if not impossible to measure the final cross- 
sectional area of the jagged fracture zone and to decide 
which projection of this area is of significance. Fur- 
thermore, the ‘“‘hypothetical ultimate’’ strain we re- 
ferred to was considered to be «& = AL/L» where Lo 
was of differential size; this value would not be the 
true strain since the latter is 


L 
=f adL/L 
Lo 


according to definition, although it is a measure of the 
same quantity. 

We have found that elongation tests are of use in the 
solution of practical forming problems when the values 
of strain determined over certain gauge lengths are 
applied to problems involving plastic deformation over 
gauge lengths of the same magnitude. Therefore, al- 
though « = AL/L» (or “engineering-strain’’) differs 
greatly from the more classic definition of strain, the 
application of the former quantity to forming problems 


€ = 
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Fic. 3. Enlarged approximately 2 diameters, 0.1 in. line grid. 










wherein the required deformations are also calculated 
as engineering strains is perfectly valid. 

The use of true strain « may be of value in translating 
stress-strain relationships in tension into stress-strain 
relationships under combined stress, such as those 
existing during the drawing process. However, as yet 
the full significance of the strain over a gauge length 
even as large as 0.01 in. has not been determined, so 
that the value of true strain or ‘‘ultimate’’ strain at the 
break is even more in question. 

If the value of true strain at the fracture should later 
prove to be of importance in flat sheet metal, it might 
well be measured through the use of a photo grid of a 
size smaller than the average grain size of the metal in 
question and substituting Zo (original grid spacing) and 
L (final grid spacing at fracture) into the true strain 
relationship. Referring to Fig. 2 of reference 2, the 
nonuniform nature of strain distribution along a speci- 
men of polycrystalline material is clearly shown in the 
outlined zone. On such a specimen a grid utilizing a 
0.04 in. gauge length would be much too coarse for an 
accurate determination of € or ¢ at the fracture; while 
a 0.01-in. grid would be reasonably accurate if we may 
assume that the grid is smaller than the average grain 
size and that the strain over a grain at the fracture zone 
is uniformly distributed. In this case, and in the case 
shown by Fig. 3, wherein a large Liider line considerably 
alters the strain-distribution pattern, it would seem 
that the photo-grid process would offer a more valid 
method of determining the true strain than the reduc- 
tion of area method. 
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A Theory for the Buckling of Thin Shells 


HSUE-SHEN TSIEN* 
California Institute of Technology 


SUMMARY 


In a series of papers written by the present author in collabora- 
tion with Th. von Karm4n and Louis G. Dunn, the effect of 
curvature of a structure on its buckling characteristics was in- 
vestigated. The purpose of these investigations was to find an 
explanation for the discrepancy between the ‘‘classic’’ theory and 
the experiments. For the case of a thin spherical shell under 
external pressure and for the case of a thin cylindric shell under 
axial compression, equilibrium states involving large deflections 
are discovered which can be maintained by loads far less than the 
so-called buckling load calculated by the classic theory of in- 
finitesimal deflections. It was felt then that since some of these 
newly found equilibrium states closely approach the observed 
phenomena the shell must “‘jump’’ suddenly trom the unbuckled 
configuration to these equilibrium states and the structures fail 
as a result of this sudden change. However, the reason why the 
shell should jump to these particular equilibrium states and not 
others was not explained. In the present paper a new principle 
involving the energy level and the geometric restraint is de- 
veloped to determine this sudden change in equilibrium states. 
By means of this new principle the buckling load of both the 
spherical shell and cylindric shell can be calculated. The agree- 
ment with experiments is good. 


INTRODUCTION 


GETHER WITH THEODORE VON KARMAN and Louis 
G. Dunn, the author" * * has suggested an explana- 
tion for the apparent discrepancy between the observed 
facts and the well-known classic theory for certain cases 
of the buckling of thin shells. The essential result of 
these investigations is the discovery of equilibrium 
states of the shell which involve a much lower applied 
load than the buckling load predicted by the classic 
theory. The load actually decreases as the deflection 
of the shell is increased. Some of the equilibrium states 
found are close to the experimentally observed ones 
both in the applied load and in the wave pattern. The 
authors then concluded that the unbuckled shell will 
suddenly jump to those equilibrium positions during the 
loading process without reaching the classic buckling 
load. However, the question why the shell should 
jump to these particular equilibrium states and not 
others remains to be answered. This is essentially the 
objection raised by K. O. Friedrichs* to the author's 
treatment of the buckling of spherical shells under ex- 
ternal pressure. In the present paper the author will 
first try to clarify this point and then apply the theory 
to a number of special cases, for which experimental 
data are available, to test the validity of the theory. 
Finally, a more precise mathematical formulation of 
Presented at the Structures Session, Tenth Annual Meeting, 
I.Ae.S., New York, January 28, 1942. 
* Research Fellow in Aeronautics. 


the buckling problem of spherical shells will be presented 
to indicate the possibilities of further investigations. 


THE CRITERIONS FOR BUCKLING 


In the author’s paper on cylindric shells under axial 
compression,’ it is stated that under average laboratory 
and actual service conditions the most probable equilib- 
rium state is the state with the lowest possible energy 
level. This statement is based upon the assumption 
that there are disturbances of sufficient magnitude so 
that the transitions from higher energy levels to lower 
energy levels are always possible. Although this as- 
sumption seems plausible, nevertheless it can only be 
verified by comparison with experiment. Since the 
available experimental data usually only give the buck- 
ling load of the specimen, the present investigation will 
be confined to buckling phenomena. 

First of all the meaning of ‘‘possible energy levels,” 
must be made clear. The two conditions for possible 
energy levels are: (1) the corresponding external forces 
and internal stresses must be in equilibrium; (2) the 
geometric restraint and loading conditions, if any, must 
be satisfied. 

The second condition depends, of course, upon the 
manner in which the specimen is loaded. To facilitate 
the mathematical treatment, two limiting cases will be 
investigated first. If the specimen is loaded in axial 
compression, then one of these limiting loading proc- 
esses is that of a rigid testing machine, and the other 
is that of ‘‘dead weight’’ loading. 

In a testing machine the parameter over which the 
operator has direct control is the distance between the 
end plates, since this distance is determined by the posi- 
tion of the loading crank. The buckling of a thin shell 
usually occurs in a fraction of a second; therefore, the 
operator really has no time to turn the loading crank 
during the jump from the unbuckled state to the 
buckled state. In other words, the distance between 
the end plates is the same at the beginning and at the 
end of the buckling process. Hence, the geometric 
restraint in this case is the end shortening of the speci- 
men. Thus, both the starting point of the jump and 
the end point of the jump must lie in a vertical line AB 
in the load vs. end-shortening diagram (Fig. 1). Since 
the end shortening of the specimen is constant, the 
potential energy of the loading force is not changed. 
The energy level of the system is given by the strain 
energy at a given value of the end shortening. Buck- 
ling will probably occur at a certain value of the end 
shortening which gives equal strain energy for both the 
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unbuckled equilibrium state and the buckled equilib- 
rium state. If these two particular equilibrium states 
are presented by points A and B, respectively, the load 
of the specimen will jump from a higher value to a lower 
value corresponding to the distance between A and B. 
The load given by A which corresponds to the un- 
buckled equilibrium state is generally recorded as the 
failing load of the specimen. 


A A 
c D 








LOAD ON SPECIMEN, P 








END SHORTENING, € 


Relation between the load and the end shortening for 
different types of loading processes. 


Fic. 1. 


In case of dead weight loading, the parameter under 
the control of operator is the amount of weight on the 
test specimen. However, he has no control of the end 
shortening of the test specimen which is determined by 
the equilibrium of the external force and the internal 
stresses. Hence, the starting point and the end point 
of the buckling process must lie on a horizontal line 
CD in the load vs. end-shortening diagram (Fig. 1). 
Because of the movement of the loading force during the 
buckling process, the potential energy of the loading 
weight must be included in calculating the energy level 
of the system. Thus, the energy level is given by the 
total potential ® of the system, 


@= S — Pe (1) 


where 5S is the strain energy, P the loading force and e 
the end shortening of the specimen. Buckling will 
probably occur at a certain value of the load P which 
gives equal total potential @ for both the unbuckled 
equilibrium state C and the buckled equilibrium state D. 
Therefore, referring to Fig. 1, the following relation 


exists 


Sc — Pec = Sp — Pep (2) 


(3) 


Sp — Sc = Pleo — €c) 


which means physically that the increase in strain 
energy of the specimen by jumping from the unbuckled 
to the buckled equilibrium states is provided by the 
decrease in potential energy of the load. 

The limiting case of a perfectly rigid testing machine 
that will apply any load to the test specimen independ- 
ent of the distance between the end plates really does 


not exist. In an actual testing machine there is always 
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a certain amount of elasticity or “give.’’ In other 
words, with the loading crank held fixed, the distanc« 
between the end plates will slightly increase if the load 
applied to the specimen is increased. Therefore, th« 
starting point and end point of jumping of a specimen 
during buckling must lie on a straight line EF in Fig. 1, 
instead of the ideal case represented by the line AB. 
The more elastic the testing machine, the less will be 
the slope of the line EF. Therefore, the effect of the 
elasticity of the testing machine is to make the buckling 
phenomena more like that which would occur with 
dead weight loading. It will be shown later that the 
failing load of a specimen depends largely on the slope 
of the line EF. Hence, the elasticity of the testing ma- 
chine has a strong influence on the failing load. 

There are certain interesting relations connected with 
the important quantities, strain energy and total poten- 
tial. Let the load P vs. end-shortening « curve of a 
structure be given by Fig. 2. Then if the structure is 
brought up to the state P:, 4, by following the equilib- 
rium curve, the strain energy stored in the structure is 
certainly equal to the external work done by the load- 
Thus, 


S(a) = [Pac (4) 
0 


which is represented by the shaded area in Fig. 2a. 
The direction followed around the contour of this area 
is as indicated and is according to the usual convention. 
From Eq. (4) it is evident that 


dS/de = P (5) 


ing force. 


In other words, the slope of S vs. ¢ curve gives the load 
on the specimen. 
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Strain energy S and total potential @ as integrals under 
the load P vs. end-shortening ¢ curve. 


Fic. 2. 






Using Eq. (1), the total potential ® is given by the 
difference of the shaded area in Fig. 2a and the rec- 
tangle OABC. In other words, it can be calculated as 
the algebraic sum of the shaded areas in Fig. 2b. Con- 
sequently, 


P: 
&(P;) = — [ver (6) 
0 


The integration is performed by following the contour 
as indicated in the figure. Eg. (6) can, of course, also 
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be obtained from Eq. (4) by partial integration. From 


Eq. (6), it is evident that 
dé/dP = —« (7) 


That is, the slope of vs. P curve gives the end short- 
ening of the specimen with a negative sign. 

Both Eqs. (5) and (6) degenerate to rather trivial 
relations in case of linear load vs. end-shortening curves, 
but they are useful in nonlinear buckling problems. 
The relations between the quantities P and ¢ and the 
two potentials S and @ are really equivalent to those 
existing between the pressure and volume and the in- 
ternal energy and the heat content in thermodynamics. 


COLUMN WITH A NONLINEAR ELASTIC LATERAL SUPPORT 


Using the results obtained in a previous investiga- 
tion,’ a preliminary calculation of the strain energy in a 
cylindric shell under axial compression shows that the 
relation between this quantity and the geometric re- 
straint—i.e., the end shortening—is rather complicated. 
It is felt that since the method of calculating the equilib- 
rium positions in this particular investigation is only 
an approximate one, the complicated relation between 
the strain energy and end shortening may not be very 
reliable. However, an exact solution of this problem is 
difficult. During the investigation of the general 
character of the buckling of shells? it was found that 
the characteristics of this type of structure can be very 
well simulated by a column (Fig. 3) supported laterally 


P 


UINITIAL POS/TION 
BUCKLED POSIT/ON 


P 


Fic. 3. Column supported laterally with a nonlinear spring. 


by a spring whose load vs. deflection curve has a de- 


creasing slope such as that shown in Fig. 4. This anal- 
ogy was first suggested by H. L. Cox. An exact solu- 
tion of the column problem can be obtained without 
any mathematical difficulty.® 

To test the validity of the theoretic analysis, the re- 
sult of calculation is compared with experimental data. 
The experimental data were obtained from a column of 
24RST of 0.090 in. thick by 0.375 in. wide and 19 in. 
long. The lateral support is a semicircular steel ring 


THIN SHELLS 





F IN LBS. 


§ IN INCHES 


Relation between the load F and the deflection 6 for a 
semicircular ring. 


Fic. 4. 


0.015 in. thick by 0.50 in. wide. This steel ring gives a 
load F vs. deflection 6 curve as indicated in Fig. 4. 
The columns were rolled to approximately a half sine 
wave. The amplitude of the initial deflection is de- 
noted by 49. The test results are shown in Fig. 5 
where Pg = r°EI//* is the Euler column load, / is the 
length, J is the moment of inertia of the column section 
and £ is Young’s modulus of the material. The dotted 
line is the theoretic curve for a perfectly straight column 
using the supporting spring characteristics shown in 






































Fic. 5. Load P vs. deflection 6 curves for columns supported 
laterally with a semicircular ring and with various amounts of 
initial deflection 5o. 
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One of the experimental curves is labeled 
However, it is practically impossible to 


Fig. 4. 
“straight.” 
obtain such a specimen; hence, the lower experimental 


values are to be expected. In Figs. 6 and 7 the maxi- 
mum loads that the columns can sustain are plotted 
against the initial deflection ratio 5)/] and against the 
deflection ratio 6// at that maximum load, respectively. 
The agreement between the theoretically calculated 
curve and the experimental values taken from Fig. 5 is 
quite satisfactory. 
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Maximum load Paz. as a function of initial deflection 


Fic. 6. 
laterally with a semicircular ring 


59 for columns supported 
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Fic. 7. Maximum load Pmar. as a function of deflection at 
the maximum load for columns supported laterally with a semi- 
circular ring. 


With this verification of the theory, a further investi- 
gation is carried out in which the relation between the 
strain energy and the end shortening of the column is the 
primary consideration. The end shortening is the 
geometric restraint for a column in a rigid testing ma- 
chine. For this calculation the load F vs. deflection 6 
relation for supporting spring is assumed to be repre- 
sented by 

F = 13.712 (1 — 208.89° + ious”) (8) 

E l l }? 
This load vs. deflection relation is similar to that shown 
in Fig. 4. The results of this calculation are presented 
in Figs. 8 and 9, where the end load ratio P/Px and the 
strain energy S are plotted against the end shortening e, 
respectively. In these figures 7 denotes the radius of 
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Fic. 8. Relation between the load P and the end shortening « 
for columns with a nonlinear lateral support and various amounts 
of initial deflection dp. 
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Fic. 9. Relation between the strain energy S and the end 


shortening « for columns with a nonlinear lateral support and 
various amounts of initial deflection 5o. 
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gyration of the column. For the case of a straight col- 
umn, the load P/Pz first follows the straight line OA 
(Fig. 8) and then the curve AB. The point A represents 
the buckled state of the column with infinitesimal ampli- 
tude 6. All the other equilibrium positions in the 
buckled state involve finite amplitudes 6. Since the 
classic theory of buckling only considers infinitesimal 
displacement, the point A can be said to represent the 
classic buckling load. 

In the strain energy diagram, Fig. 9, the correspond- 
ing points are denoted by the same letters as in Fig. 8. 
Thus, the straight equilibrium states of the column are 
represented by the parabola OA, and the buckled equilib- 
rium states by the curve ABCD. It is interesting to 
notice that the curve AB closely follows the parabola 
OA—.e., the strain energy is only slightly higher than 
that of straight column with the same end shortening. 
The curve turns sharply at B and meets the parabola 
at C. The branch CD has less strain energy for the 
same end shortening than either the straight column or 
the buckled state represented by AB. Therefore, ac- 
cording to the criterion of lowest strain energy, the 
actual behavior of the column in a testing machine 
when the end plates are gradually brought together will 
be, first, the straight unbuckled shape up to C. AtC 
there will be a sudden change in the shape of the column 
to that which involves a finite amplitude 6 and a jump 
in the magnitude of the load from P/Pz = 2.40 to 
P/P: = 1.09. Since, according to Eq. (5), the slope 
of the S vs. e curve is a measure of P/Pz, this jump in 
the magnitude of P/P: is indicated in Fig. 8 by a sud- 
den change in the slope of the two branches of the curve. 
From then on the column will behave according to the 
curve CD with increasing end shortening. Therefore, 
under ordinary conditions the classic buckling state 
does not appear in the buckling process at all, and the 
recorded failing load of the column will be P = 2.40P: 
instead of the classic value of P = 3.61Px correspond- 
ing to the point A. 

The initial deflection 59 of the column tends first to 
decrease the failing load. Finally, the sudden jump 
in the value of P disappears completely when, at any 
given value of end shortening e, there is only one value 
of the load P. For this particular problem the small 
initial deflection of 6) = 0.17i is already big enough to 
eliminate the jump in the value of P (Figs. 8 and 9). 

To investigate the behavior of the columns when 
loaded by weights, the total potential has to be plotted 
against the load P. This is done in Fig. 10, using the 
same characteristics of the lateral support as given by 
Eq. (8). For the straight column, 6)/7i = 0, the un- 
buckled equilibrium states are represented by the 
parabola OCoAo. The classic buckling load is indicated 
by the point Ap. The curve ApBoCo represents the 
buckled equilibrium states. Using the buckling cri- 
terion developed in the previous sections, the column 
will maintain its unbuckled straight shape up to the 
point Cy with P/P, = 1.196, while the loading is gradu- 
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Fic. 10. Relation between the total potential and the load 
P for columns with a nonlinear lateral support and various 
amounts of initial deflection 59 


ally increased. At this value of P/P: the column will 
suddenly jump from the unbuckled equilibrium state 
to the buckled equilibrium state. At the same time the 
value of the end shortening will change from ¢/] = 
1.199?(z/1)? to e/l = 5.60r*(i/1)*. According to Eq. 
(7) this jump in the magnitude of the end shortening is 
signified in Fig. 10 by the sudden change in the slope of 
the two branches of the curve. Here, again, the failing 
load will be P = 1.196Pz¢ instead of the classic buckling 
load of P = 3.61P:x. 

The initial deflection 49 of the column again tends first 
to decrease the failing load. For 69/71 = 0.1 the failing 
load is P = 1.088Pz; for 69/7i = 0.5 the failing load 
is P = 0.778Pz. Finally, when at any given value of P 
there is only one equilibrium value of ¢, the jump in the 
equilibrium states of the column during the loading 
process will disappear completely. 

Summarizing the results, the failing load of a column 
with a nonlinear elastic support as given by Eq. (8) is 
given in Table 1. It is seen that under ordinary condi- 
tions the so-called classic buckling load really does not 
necessarily indicate the load carrying ability of the 


column. 


TABLE 1 
P/Px —P/P,» of New Theory 
of Classic Testing Dead 
bo/mi Theory machine weight 
0 3.61 2.40 1.196 
0.1 1.79 1.79 1.088 
0.5 0.778 


| & 


0.88 


THIN UNIFORM CYLINDRIC SHELL UNDER AXIAL 
COMPRESSION 


With the buckling criterions so clearly demonstrated 
by the failure of a column with nonlinear elastic sup- 
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port, the theory can now be applied with confidence to 
thin cylindric shells under axial compression. Using 
the results already obtained by the author in a previous 
paper,® the strain energy S of the shell and the total 
potential can be calculated. The results are shown 
in Figs. 11 to 14. In these figures R is the radius, ¢ is 
the thickness and ¢ is the unit end shortening of the 
shell. Furthermore, s denotes the strain energy per 
unit area of the shell surface, and ¢ is the total potential 
per unit area of the shell surface. The number of 
waves m in the circumference of the shell is determined 
by the parameter 7 which is equal to n*(t/R). The 
ratio of the wave length in the circumferential direction 
to that in the axial direction is defined as u. Thus, for 
uw = | the waves are square; for u = 0.5, the waves are 
rectangular with long side parallel to the axis of the 
cylindric shell. The value of Poisson’s ratio is taken to 
be 0.3. 
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Fic. 11. Relation between the unjt strain energy s and the 
unit end shortening e for cylindric shells under axial compression 
with the aspect ratio u of the waves equal to 1.0. 


To investigate the behavior of the shell under a rigid 
testing machine, the strain energy parameter s(R/t)?/Et 
is plotted against the end-shortening parameter «R/t as 
indicated in Figs. 11 and 12. According to the criterion 
developed in the previous sections, buckling will prob- 
ably occur when the strain energy curve for the buckled 
equilibrium positions first crosses the curve for the un- 
buckled equilibrium positions. Thus, by referring to 
Fig. 11 for » = 1.0, the lowest probable buckling load 
is indicated by the point A. Then the corresponding 
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Fic. 12. Relation between the unit strain energy s and the unit 
end shortening ¢« for cylindric shells under axial compression 
with the aspect ratio » of the waves equal to 0.5. 
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Fic. 13. Relation between the unit total potential g and the 
compression stress o for cylindric shells under axial compression 
with the aspect ratio u of the waves equal to 1.0. 
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buckling load can be calculated as oR/Et = 0.460 where 
o is the average axial compression stress. The value of 
nis 0.400. Thus if R/t = 1000,” = 20. At buckling, 
the stress will jump from ¢ = OA60Et/R to o = 
0.239Et/R. With uw = 0.5, Fig. 12 shows that the low- 
est probable buckling stress is around 


o = 0.370Et/R (9) 


The corresponding value of 7 is 0.289. Thus if R/t = 
1000, » = 17. At buckling, the stress will jump from 
the value given by Eq. (9) to o = 0.2000E#/R. 

The failing load of the cylindric shell under dead 
weight loading can be determined from Figs. 13 and 14, 
where the total potential parameter y(R/t)*/Et is 
plotted against the stress parameter ¢R/Et for » = 1.0 
and uw = 0.5, respectively. Using the criterion that 
buckling will probably occur for dead weight loading 
when the total potential curve for the buckled equilib- 
rium state crosses the curve for the unbuckled equilib- 
rium state. For uw = 1.0, this happens at point A 
(Fig. 13), which corresponds to ¢ = 0.298Et/R. Again, 
since there is a sudden change in slope of the two 
branches of the curve which the shell will follow, the 
buckling phenomena is accompanied by a sudden in- 
crease in end shortening. By Fig. 14, for u = 0.5, the 


lowest probable buckling stress is around 
o = 0.238Et/R (10) 


Again, there will be a sudden change in the end shorten- 
ing at buckling. 
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Fic. 14. Relation between the unit total potential y and the 
compression stress o for cylindric shells under axial compression 
With the aspect ratio u of the waves equal to 0.5. 
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Both Eqs. (9) and (10) give a buckling stress much 
less than the classic buckling stress of ¢ = 0.606Et/R. 
As shown in the previous sections, the elasticity of the 
testing machine tends to lower the buckling stress 
toward that which would be obtained by dead weight 
loading. However, even the dead weight failing load 
given by Eq. (10) is higher than the average of the 
reported experimental data which give o = 0.15Et/R.? 
One of the reasons for this failure of quantitative agree- 
ment between the theory and experiment is, of course, 
the approximate character of the method used in refer- 
ence 3 which is the basis of the present calculation. 
Furthermore, the aspect ratio wu of the wave is chosen 
rather arbitrarily. There may well be other values of u 
which will give lower values of the buckling stress. 
These questions can only be clarified by investigations 
to obtain more accurate solutions of the problem. 


SPHERICAL SHELLS UNDER EXTERNAL PRESSURE 


The buckling of thin uniform spherical shells under 
external pressure was investigated by the author in 
collaboration with Th. von Karman! and by K. O. 
Friedrichs.‘ In the author’s treatment there are sev- 
eral simplifying assumptions that were shown by 
Friedrichs to be unnecessary. Friedrichs, however, 
was unable to give a definite answer to the problem. 
Therefore, this problem will be treated now in the light 
of the present theory. 

It was shown by Friedrichs that if the ratio of thick- 
ness to the radius of the shell and the angular extension 
of the buckle is small and if the additional radial and 
circumferential stresses due to buckling are zero at the 
boundary of the buckle then the difference between the 
sum of the additional strain energies of extension and 
bending and the virtual work of the external pressure 
force acting on the surface of the shell can be calculated 
by the following integral: 


(+/R)2 B/a. \2 
I= :_¥ R)* = (=) akda + 


2 12(1 — v?) 0 \da 


8 8 , “8 
cf | f (: — vada ada = = ri aida (11) 
2/0 @ 2 2EJ/o 


Here, ¢ is the thickness, R is the radius of the spherical 

shell, 8 is the angular extension of buckle as shown in 

Fig. 15 and o is the stress in the unbuckled shell under 

an external pressure , and is equal to 

o = PR/2 (12) 

If, due to buckling, the inward radial displacement of 
the shell element is w, the function x(a) is given by 

x(a) = —(1/Ra)(dw/da) (13) 

To determine the equilibrium states the energy 

method can be used. For this purpose the radial dis- 

placement w is assumed to be 


w/R = f[1 — (a/B)*]* (14) 
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Fic. 15. Symmetric buckling of a spherical shell segment under 


external pressure. 


where f is the ratio of the maximum radial displacement 
at the center of the buckle to the radius of the shell. 
This particular form of radial displacement has a zero 
slope at the boundary a = 8. However, the curvature 
is not zero at this boundary. Therefore, it requires a 
distribution of bending moment at the boundary to 
keep the buckled shell segment in equilibrium. Substi- 
tuting Eq. (14) into Eq. (11), the integral J can be 


written as 
4 (¢/R2 ff, B/1., 2f, &ft\ lao 
Sea > oe.¢ a Be ng i i ly 
9 (1 — »v?)B? . 2 os f 9 B? . 63 a) 3) 
(15) 


To find the value of f at which the shell can be in 
equilibrium with a given value of o, the derivative of J 





with respect to f has to be set equal to zero. Thus 
2 ¢ 2 2 
peta &+a5) 
E 3(1 — vp? 20 26? 21 6 


Eq. (16) can be rewritten in a simpler form by using the 
following parameter: 


t= BR/t, § = fR/t (17) 
Then, 
Pee ot ee ee 1. 8# 
~ a GE Se ees ey, 
HE 30-202 2° * a7 (18) 


When the amplitude of the buckle is very small, i.e., 
& — 0, the minimum value of the buckling stress is given 
by 


o = 0.9376Et/R (19) 


with vy = 0.3. This value is higher than the classic 


value 


« = 0.606Et/R (20) 


due to the special assumed deflection form given by 
Eq. (14). However, for the present calculation, the 
values of oR/Et for finite values of the amplitude ratio & 
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are more interesting. Using Eq. (18), the value of 
oR/Et can be calculated as a function of £ for various 
fixed values of ¢ which determines the size of the 
buckle. The results are shown in Fig. 16. 





ef 


Fic. 16. Relation between the compression stress o and the 
maximum deflection 6 for a spherical shell segment under ex- 
ternal pressure. 


The envelope of curves shown in Fig. 16 can be ob- 
tained by differentiating Eq. (18) with respect to ¢ and 
then setting the resultant expression equal to zero. 
Thus, 

SF. 


21 ¢? 


4 1 1 3 
cw Senta geile (21) 

3(1 — »)¢2 20 
Combining Eqs. (21) with (18), the equation for the 
envelope curve is obtained as 


R 1/16 1 32 1 
(<*) = 1 ge — + = --& (22) 
Et) mun 2°5(1—») 35 2 


For large values of ¢, the envelope value of ¢R/Et will 
become negative and decrease without limit. This is 
different from that which the author found in his previ- 
ous investigation! where a minimum value of oR/Et 
exists. This minimum value was due to the particular 
displacement of the shell element used in that investiga- 
tion. The negative values of ¢R/ Et at large deflections, 
however, are not the failing stress of the shell, as will 
be shown presently. 

When a segment of a spherical shell is loaded by ex- 
ternal pressure, the experiments show that in general 
only one single small buckle appears, whereas the rest 
of the shell remains in approximately spherical shape. 
Therefore, in calculating the strain energy and total 
potential of the shell, both the buckled portion and the 
unbuckled portion have to be considered. In the pres- 
ent calculation the radial displacement of the shell ele- 
ment is assumed to be composed of two parts: the uni- 
form inward displacement equal to (1 — v)Ro/E due 
to the external pressure and an additional displacement 
given by Eq. (14) within the boundary of the buckled 
portion. Furthermore, the additional extensional 
stresses due to buckling are assumed to be zero at the 
boundary of the buckle. These assumptions result in 
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the following inaccuracies: (1) There is a discontinuity 
of curvature of the shell surface at the boundary of the 
buckle; (2) there is a discontinuity of the circumferen- 
tial displacement of the shell elements at the boundary 
of the buckle; (3) the effective edge support at the 
boundary of the shell segment is neglected. 

Since in all the calculations the ratio of the thickness 
to the radius of the shell is taken to be very small, it is 
probable that the errors introduced by the above- 
mentioned inaccuracies will not be appreciable. 

With these assumptions, the strain energy s per unit 
area of the shell surface can be calculated by using Eq. 
(14) and expressions similar to Eq. (11). The result is 


sR? <4) a E oR 
a hoe 1S > 
(1 — v)#E (A R(i— v) L3 Et 


4 1 ] 1oR 1 ' “4 
( Sa ie, xt) —* + 4 (233) 
91 —»)c? 20 3£Ett/” 9 63etl' 





where \ is the ratio of area of a hemisphere with radius 
R to the area of spherical segment. Therefore, if 6 is 
the semiangular extension of the shell segment, A can be 
written as 

(24) 


4 = 1/(1 — cos 6) 


If the shell is tested by means of external fluid pres- 
sure and the fluid is incompressible and enclosed in a 
chamber, the geometric restraint is the volume change 
of the shell due to its deflection. This volume change 
is caused by the uniform inward displacement due to the 
external pressure and the additional displacement given 
by Eq. (14) within the boundary of the buckled region. 
The change of volume v per unit area of the shell surface 
can be shown to be 


v _ oR 


3 , L/RIEE 
(l—»v)t Et 


6(1 — »v) 


(25) 


This type of closed chamber fluid pressure test is equiva- 
lent to the so-called rigid testing machine test previ- 
ously treated. From the similarity to the previously 
treated cases, the probable failing load of the shell is 
determined by the condition that the unbuckled and 
buckled equilibrium states have the same strain energy 
and same volume change. That is, if the quantities 
corresponding to the starting point of buckling are de- 
noted by the subscript 1 and those corresponding to the 
end point of buckling by subscript 2, the conditions are 


i =& =0 (26) 
OR _ oR, yt _ tbs (27) 
Et Et R 6(1 — v) 
ak (28) t Sake {i o2R 
=(—)} +\,_—- {- = + 
( t R(i-—yv)'3 Et 
4 1 1 1 oR 1 te 4 &?] 
—— +— --= |g -2= 4+ —= 
91 — py?) f2? 20 3 Et fe Of 63 $2" 


(28) 
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Eliminating o, from Eqs. (27) and (28), the following 
relation is obtained. 


a a ye 
eg = lle & + 
63 ¢2" Ifo 
a l 1 2 ¢ 
ie +4 ~ 1 s l ax A(t R)& at 
11=)> ff DW ££ 36(1 — pv) 


(29) 


However, since at the end point of the buckling process 
the shell is at equilibrium, o2, f2 and & must also satisfy 
the equilibrium condition as given by Eq. (18). Thus, 


4 ‘ » 
3 3 J ~~ 


8 l 
3(1 — v?) fe 20 Et 


Ste l 
2” aal 
(30) 
Eqs. (29) and (30) give o2.R/Et and & as functions of 
¢2 for a given value of \4/R. Then from Eq. (27) the 
expressions for o,R/Et and o2R/Et are, 
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Fig. 17 indicates the calculated results for ¢/R = 
0.02101. It is seen that as f2 increases the value of 
o1R/Et first decreases and then increases again. The 
curves terminate at a value of f2 which causes the ex- 
pression under the radical in Eqs. (31), (32) and (33) 
to be equal to zero, because any further increase in {> 
will give imaginary values to the radical. 

Since for a given shell segment only the value of \//R 
is fixed, the size of the buckle or the parameter {2 is free 
to vary, except that the size of the buckle cannot be 
larger than the size of the segment. Assuming for the 
moment that the last condition is not violated, then the 
actual probable failing load is determined by the mini- 
mum value of o,R/Et. Denoting this minimum value 
by o1*R/Et and the corresponding o2R/Et by o2*R/Et, 
then from Fig. 17, for 44/R = 0.02101, o«*R/Et = 
0.320 and o2*R/Et = 0.125. 

By carrying out a series of computations for different 
values of \t/R, o1;*R/Et and o2*R/Et can be plotted 
against the parameter R/t\. This is shown in Fig. 18. 
The upper curve which represents o,*R/Et then gives 
the failing load of the shell segment under the previ- 
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Fic. 17. Initial buckling stress o; and final buckling stress o2 


of a spherical shell segment under external pressure as a function 
of the angular extension 8 of the buckle. 
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Fic. 18. Initial buckling stress o,* and final buckling stress o2* 
of spherical shell segments under external pressure. 
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Fic. 19. 





Angular extension £,* of the buckle for spherical shell 
segments under external pressure. 
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ously stated conditions. This figure shows that the 
failing stress parameter o,*R/Et is no longer a con- 
stant as predicted by the classic theory but decreases 


‘slowly with increasing values of R/t if the solid angle of 


the segment—.e., \X—is fixed. In fact, at a fixed value 


of \ it is found that 


o1*/E ~ (t/R)'-* (34) 


for 20 < R/dt < 100. Furthermore, as shown in Fig. 
19, the value of {2* at which o,*R/Et occurs also de- 
creases with increasing R/t. It is found that with a 
fixed value of \ within the range stated {2* varies ap- 
proximately as (R/t)°-*8. Then, according to Eq. (17), 
62* varies approximately as (R/t)~°-*°. However, if 
lis the diameter or wave length of the buckle, //2R = 
Bo*. Then 


1/t ~ (t/R)~%-*4 (35) 


The exponents of ¢/R for a buckling stress ¢/E and a 
wave length //R given by the classic theory are, of 
course, 1 and —!/, respectively. 

There is no experiment known to the author which is 
performed under the geometric restraint that the en- 
closed volume of the shell during the buckling must be 
constant. However, a number of experiments have 
been conducted at the California Institute of Tech- 
nology on clamped spherical shell segments of 6 = 17 
deg. 45 min. subjected to water pressure.t The water 
is slowly admitted to a closed chamber whose top cover 
is the shell segment. The water is immediately turned 
off as soon as the shell starts to buckle, so that the vol- 
ume of water in the chamber is nearly constant during 
the buckling process. However, the pressure of the 
water is measured by a mercury manometer connected 
to the chamber by a piece of rubber tubing which is 
filled with air at the beginning of the test. Both the 
drop in height of the mercury column and the expansion 
of air trapped in the rubber tubing will effectively in- 
crease the volume of the chamber. Therefore, these 
tests are somewhat similar to that of a column loaded 
in a partially elastic testing machine. The results of 
these tests are shown in Fig. 20 where the values of 
oR/Et for the beginning of buckling and the end point 
of buckling are plotted against (R/t). It is seen that 
both values tend to drop off slightly as R/t increases, 
which is a feature also shown by the theoretic calcula- 
tion for constant volume restraint. The dotted curves 
represent o,*R/Et and o2*R/Et previously calculated. 
Although the test points do follow the curves very well, 
the agreement can only be regarded as qualitative be- 
cause of the difference in the test conditions of these two 
cases. 

If the shell segment is loaded by means of fluid pres- 
sure and the fluid is drawn from a large reservoir with a 
free surface, then the small change of volume during the 


{+ The experiments were carried out by Mr. H. B. Crockett 
under Dr. L. G. Dunn’s supervision. The author is deeply 
grateful to both of them. 
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Fic. 20. Initial buckling stress and final buckling stress of 
spherical shell segments loaded with water pressure (experiments). 


buckling process will not appreciably change the level 
of the fluid in the reservoir. Under such test conditions, 
the pressure acting on the shell at the starting point of 
buckling is practically the same as that at the end point 
of buckling. Let the volume change be AV: then the 
corresponding work done by the fluid is given by pAV. 
This amount of work done must be equal to the in- 
crease in strain energy of the shell. This will give a 
similar relation as Eq. (3). Therefore, if the total 
potential @ of the system is defined as 


&’= S$ — pv (36) 
the buckling will occur when, at equal values of oR/Et, 
the total potential of the unbuckled equilibrium state 
is the same as that of the buckled equilibrium state. 
By means of Eqs. (36), (23) and (25), the total potential 
¢g per unit area of the shell surface can be expressed 
as 


nR2 2 = £2 
oR? iid (=) 6 »/ Le x 
(1 — v)#E Et R (1 — v) 


4 1 1 1oR 1 & 4 = 
{| + is re |- = t (37) 
(Loi —»v)¢ 20 3&¢ Be 
Let o3;3R/Et be the buckling stress parameter for this 
type of loading; then the buckling criterion stated above 
requires that 


— (EY = — (22) ae af Saks" x 
Et Et R (1 am v) 
‘| Bie ha des ce BE ‘] ~ 4 4 &'L (gg) 
(Lol — »*) fs? 20 3 Ethst 9% 63 gS 
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e ki 1 & 
ted a 3 
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1 o3R l 


- = 0 (39 
3 Et | 


| 4 1 ee 
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On the other hand, since the buckled shell must be in 
equilibrium, Eq. (18) must also be satisfied. Therefore, 
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From Egs. (39) and (40) o;R/Et and &; can be obtained 


as functions of {;3. Then, 
R 1 
a: -< O E 7 td (41) 
Et 240 3(1 — »?) ts 
and 
& = (7/8)é3 (42) 


If the shell segment is large enough, the actual buck- 
ling load will be determined by the minimum of o;R/E£t. 
Denoting quantities corresponding to this minimum by 
an asterisk, calculations show that 


o3*R/Et = 1/3V5(1 — v*) (43) 
f3* = 8V5/(1 — y2) (44) 

and 
§* = 7V5/(1 — v?) (45) 


It should be noted that under this type of testing the 
buckling stress parameter o;*R/Et is again a constant 
However, the numerical value 
of this constant is only about one-fourth of the classic 
value. 

The experiment on a hemispherical shell made by 
E. E. Sechler and W. Bollay was performed by immers- 
ing the shell in a mercury bath. Hence, the test condi- 
tions approximate the constant pressure type to which 
Eqs. (43), (44) and (45) applies. In Table 2 the values 
predicted by Eqs. (43), (44) and (45) are compared with 
those deduced from the test. The agreement between 
theory and experiment is truly remarkable. 


as in the classic theory. 


TABLE 2 


Experiment Theory 
o3R/Et 0.154 0.1561 
B;* 8 deg. 8.3 deg. 
f;* = 8*/t 12.5 16.4 


CONCLUDING REMARKS 


In the previous sections a new principle for determin- 
ing the buckling load of curved shell under ordinary test 
and service conditions has been developed. This prin- 
ciple of lowest energy level is verified by comparing the 
experimental data with the theoretic predictions. 
However, in view of the prerequisite that arbitrary dis- 
turbances of finite magnitude have to exist, the buckling 
load determined by this principle may be called the 
“lower buckling load.’’ The classic buckling load that 
assumes only the existence of disturbances of infinitesi- 
mal magnitude may be called the “upper buckling 
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load.’”’ Of course, by extreme care in avoiding all dis- 
turbances during the test, the upper buckling load can 
be approached. The lower buckling load, however, has 
to be used as the correct basis for design. 

It may well be mentioned here that the possibility 
of the existence of the lower buckling load is connected 
with the elastic characteristics of the structure—i.e., 
the load on the structure must decrease as the deflection 
is increased. If this is not the case, the lower buckling 
load does not occur, and the buckling load of the struc- 
ture can be closely calculated by the classic theory, as 
shown by the well-known case of a flat plate under edge 
compression. The author believes that when a cylin- 
dric shell is loaded by external lateral pressure, the load 
rises with the increase in the deflection of the shell sur- 
face after buckling. Therefore, for this case, also, the 
buckling load can be closely predicted by means of the 
classic theory of thin shells.’ 

To indicate the possible future developments, the 
problem of buckling of a spherical shell of constant 
thickness under uniform dead weight pressures will be 
formulated more precisely. If the radial and circum- 
ferential displacements of the median surface of the 
shell from its spherical equilibrium positions are de- 
noted by w and u, respectively, the strains and change 
in curvatures can be calculated as‘ 


¢, = Gu/da) — w , 1 [(dw/de) + u}? 
R 2 =R 


(u cota — w)/R 
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a is the angle between the radius vector and the axis of 
symmetry. The difference in the total potential  — 4, 
for the buckled and unbuckled state with same external 
pressure can then be expressed as 


[ («2 + «7 + 2ve.) 
/S0 
iiaiie a t? fi 
sin —____— (x7 7+ 2vxx, (47 
ce 12(1 — »?) , Ke Xt H 20Kex my 
sin ada — ead ray tu) R? -— 
E Jo 2 R? 
t ws +(e )- 
tee = — — pote 
3R? R R 
w* 


zs —8(1 — ed sin ada 


b’-% 1 
wR°Et (1 — vp?) 





All equilibrium positions can be determined by the 
condition that the variations of the integral with respect 


to wand u must be zero. Or 
6[(®@ — %/xR°Et)|,, = 0 (48) 
5[(@ — &/rR*Et)|,, = 0 


The upper buckling load is determined by the additional 
condition that 


uO (49) 


WwW = 


The lower buckling load is determined by the additional 
condition that 


(@ — &)/rR°Et = 0 (50) 


The problem of determining the lower buckling load 
thus becomes a special case of the well-known isometric 
problem in variational calculus. 
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Tricycle Landing Gear Design—Part 1° 
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ABSTRACT 


6 hen PROBLEMS OF THE tricycle landing gear are dis- 
cussed in the light of available information with 
the object of providing criteria to assist the designer. 
The general geometric arrangement involving the deter- 
mination of wheel base, tread, and center of gravity 
location is first considered. It is shown that the nose 
wheel should be located as far forward of the center of 
gravity as possible, and that the fore-and-aft location 
of the rear wheels is limited to a narrow range by condi- 
tions of balance and longitudinal stability. The rela- 
tionship between tread, wheel base, and the resistance 
to turning over is found, and the effects of tread and 
fore-and-aft location of the rear wheels on the direc- 
tional stability and ground maneuverability are dis- 
cussed. 

The problems related specifically to the design of the 
nose wheel are next examined, and a basis for the selec- 
tion of nose wheel and tire size is given. The fundamen- 
tal causes of shimmy are reviewed, including the 


effects of trail, caster angle, wheel offset, tire type, and 
moment of inertia of castering parts. Shimmy elimina- 
tion is discussed with special reference to elimination 
by damping and including methods of calculating the 
amount of damping ‘necessary to prevent shimmy. 
The construction of fluid dampeners is described, and 
their damping characteristics are compared to the 
simpler mathematically defined types of damping. An 
empirical relationship between the volume of a vane- 
type hydraulic dampener necessary to prevent shimmy 
and the static nose wheel load is given. 

Ground handling and maneuverability, construction 
and installation difficulties, castering locks and stops, 
and steering devices are briefly discussed. Conclusions 
as to the best arrangements must be based on the char- 
acteristics of the particular design. As a partial check 
on the validity of the criteria developed, the results of 
applying them to some successful tricycle gear airplanes 
are summarized. 





INTRODUCTION 


The susceptibility of the conventional tail wheel 
landing gear to nosing over, its inherent tendency to 
ground loop (particularly in a side wind), and the 
necessity of making a full-stall landing to prevent 
bouncing have made it less and less satisfactory as 
landing speeds increased. The tricycle landing gear 
will eliminate these undesirable characteristics if 
properly designed, but it may introduce other objection- 
able features peculiar to its type. Some of these possible 
difficulties are porpoising, longitudinal instability with 
the nose wheel off the ground, shimmy of the nose 
wheel, and failure of the nose wheel strut from excessive 
dynamic loads. In addition, it should be remembered 
that the very fact that the faults of the conventional 
gear are eliminated will result in increased loads on the 
structure. Removal of the possibility of nosing over 
will lead to faster taxiing and use of full power on muddy 
fields, elimination of the ground-loop tendency will re- 
sult in more cross-wind landings, and the lack of any 
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tendency to bounce will mean more landings at high 
sinking speeds. 

Many investigations, both theoretic and experi- 
mental, have been made of tricycle gear characteristics. 
Most of these, however, are either of.a rather general 
nature or do not present their results in a manner that 
may be easily used in design work. Frequently, the 
conclusions reached are erroneous as a result of simpli- 
fying assumptions that ignore important factors. Even 
design criteria of the various authorities are sometimes 
unreliable since they are usually only an adaptation of 
the conventional gear criteria to the geometry of the 
tricycle type or have not been revised in the light of 
practical experience. 

This paper is an attempt to provide reliable informa- 
tion on the tricycle gear obtained from both the litera- 
ture on the subject and independent investigations. 
Considerable attention has been given to presenting the 
information in a manner to permit the quick solution 
of the various design problems. Most of the results 
are presented graphically in a series of charts which 
may be applied directly in the majority of cases. The 
paper is separated into two major divisions. Part 1 
deals with the general geometry of the gear, such as the 
wheel base, tread, and location relative to the airplane. 
Part 2 is concerned with ,the nose wheel design and 
considers the problems of nose wheel size, trail, caster 
angle, and shimmy dampeners. 
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(GENERAL (GEOMETRY 


WHEEL BASE 
General 


The effect of wheel base on the characteristics of the 
landing gear cannot be determined without considera- 
tion of its locations relative to the center of gravity and 
mean aerodynamic center of the airplane. The analysis 
becomes much simpler if instead of considering the 
effect of wheel base as a whole it is broken down into 
the effect of the horizontal distances of the front and 
rear wheels from the c.g. and the distance of the rear 
wheels aft of the aerodynamic center. 


Nose Wheel Location 


The principal item affected by the distance of the 
nose wheel forward of the c.g. is the nose wheel load. 
Other characteristics related to this dimension are the 
stability against overturning and the resistance to 
porpoising. 

In attempting to find the landing loads on the nose 
wheel, some consideration must be given to the possible 
landing attitudes of the airplane. The three principal 
landing modes for a tricycle gear type may be defined 
as three-wheel landing, rear wheels first landing, and 
nose wheel first landing. In the three-wheel landing, 
the load distribution between the front and rear wheels 
is usually assumed to be the same as the static distribu- 
tion. Actually, the shock absorber characteristics may 
result in a distribution considerably different from the 
static. By the conventional assumption, however, the 
nose wheel load in a three-wheel landing is 


Py = nW(l — 1,)/l (1) 


where /7 is the wheel base; /;, the horizontal arm from 
the c. g. to the nose wheel; m, the acceleration; and W, 
the airplane weight. 

In a rear wheels first landing, the proximity of the 
rear wheels to the c.g. will result in most of the energy 
being absorbed there, after which the airplane will 
drop over on the nose wheel. The nose wheel load from 
this condition may therefore be expected to be less 
than in the three-wheel landing. 

In a nose wheel first landing, the nose wheel striking 
the ground tends to rotate the airplane. The load on 
the nose wheel thus depends on the pitching moment 
of inertia of the airplane. This load can be determined 
in terms of the acceleration at the nose wheel. The 
acceleration at the nose wheel may be found, in turn, 
from the shock absorber travel and airplane sinking 
speed. 

Let 


horizontal distauce from the nose wheel to 
the c.g. 
airplane weight 


I; = 


VW = 





W/g = airplane mass 

g = acceleration of gravity 

k = airplane pitching radius of gyration 

a = vertical acceleration at airplane c.g. 

dy = vertical acceleration of the airplane at the 
nose wheel 

w = angular acceleration of airplane 


If the airplane strikes the ground on the nose wheel 
with a force, F, the torque about the c.g. is Fl; and the 
angular acceleration is the torque divided by the pitch- 
ing moment of inertia or 


w = Fl,/k?(W/g) (2) 


The acceleration at the c.g. produces the force F, 
or 
F = Ma = (W/g)a (3) 


The acceleration of a point on the airplane above 
the nose wheel is equal to the acceleration at the c.g. 
plus the effect of the angular acceleration. The angular 
acceleration has a vertical component at the nose wheel 
of wl;. The total vertical acceleration of the airplane at 
the nose wheel is therefore 


dy =~at+oal, (4) 
Substituting Eq. (1) in Eq. (2) 
ay = a+ [FI,*/k?(W/g)] 
Replacing F by Eq. (3) 


dy = a+ (W/g)al,?/(W/g)k? 
or 

dy = a(l + (1,2/k?)] 
Therefore 


a = ay[k?/(R? + 1,”) | (5) 


The acceleration at the airplane c.g. is thus always 
less than at the nose wheel except when /, becomes zero 
or what amounts to the same thing if the nose wheel 
were put at the c.g. In this special case, of course, the 
acceleration at the nose wheel becomes the same as the 
acceleration at the c.g. 

The load on the nose wheel, F, can now be found by 
use of Eq. (3). 


F = (W/g)a = (W/g)an[k?/(R? + 1,°)] (6) 


If the acceleration is given in terms of gravity, the 
term dy/g can be replaced by ny 

F = Wny[k?/(k? + 1,7)] (7) 

the 


The acceleration, my, acting on the airplane at 
nose wheel is determined by the vertical velocity at 
which the airplane strikes the ground and the char- 
acteristics of the nose wheel shock absorber. The value 
of my can be determined by a drop test of the nose wheel 
assembly from a height that will give the desired sinking 
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speed at the time of impact with an equivalent weight 
on the nose wheel of 


w= W(k/(k? + 1%] (8) 


The load factor determined from this drop test is my and 
must not exceed the design load factor for the nose 
wheel and supporting structure. If the nose wheel 
travel and the sinking speed are the same as that of 
the main gear, the load factor my will be approximately 
the same as that of the main gear in the three-wheel 
landing condition. Should the nose wheel travel be 
different from that of the rear wheels, the nose wheel 
load factor would be affected approximately in the 
inverse ratio of the wheel strokes. As the nose wheel 
travel is usually the same as the main wheels and as 
the same sinking speed is usually assumed for all condi- 
tions, the load on the nose wheel in a nose wheel first 
landing may be taken as 


Py = nW{[k?/(k? + h2)] (9) 


where m is the load factor for a three-wheel landing. 

From Eq. (9) it is apparent that, unlike the three- 
wheel landing where the nose wheel load is expressed 
by Eq. (1), the total wheel base has no effect on the 
front wheel reaction. The term »W appears in both 
equations, but the percentage of »W acting on the 
nose wheel is determined by entirely different factors. 
In the three-wheel landing, the distribution, being done 
by simple statics, is easily understood. The meaning 
of the expression k?/(k? + 1,7) in Eq. (5) through Eq. 
(9) is not so easily recognized. A simple analogy will 
help explain its significance. 

The terms in Eq. (1) are unaffected by the distribu- 
tion of the items comprising the weight, W, with re- 
spect to the c.g. so long as the c.g. remains at the same 
point. Thus, if the engine were taken out of the nose 
of an airplane and put near the c.g. position and the 
radiators and some other equipment were taken from 
the tail and moved far enough forward to balance the 
airplane at the original c.g. position, the nose wheel 
load as determined by Eq. (1) for a three-wheel landing 
would remain the same. The load for a one-wheel land- 
ing would be markedly reduced, however, because the 
shift of the masses making up the airplane would have 
reduced k*. Thus, although the c.g. and the weight 
remained the same, the nose wheel load was reduced. 
With the engine in the nose immediately above the 
nose wheel, the engine inertia had to be absorbed by 
the nose wheel. Moving the engine back in the air- 
plane left the engine practically free to continue fall- 
ing when the nose wheel hit and thus largely relieved 
the nose wheel of decelerating the engine. It is possible 
to represent the airplane by an arrangement of two 
concentrated masses that will produce the same loads 
on the landing gear. In this equivalent airplane (Fig. 1) 
the weight is divided so that a concentrated weight 
equal to Eq. (8) is located directly over the nose wheel. 
The remainder of the weight (W — w) is concentrated 
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Fic. 1. Equivalent weight distribution. 


at a distance k*//, aft of the c.g. position. The c.g. 
therefore remains in the correct position, as, taking 
moments about the c.g. of the two concentrated 


weights, 


Moment = w(i* ) ; w( a) wie, 
Rk? + 1,3 k? + 1?/ |, 


w (FA =) =" 
k + I? 


Therefore, the wheel loads in a three-wheel landing are 
given by Eq. (1). 

In a nose wheel first landing, the nose wheel has to 
stop only the weight, w, immediately over it. The re- 
mainder of the weight, (W — w), is free to continue 
falling until the rear wheels strike the ground. If the 
airplane is taxiing and the nose wheel strikes a bump, 
the nose wheel must lift the weight, w, over the bump. 
The distribution of wheel load by statics would hold 
only if the rear wheels struck a bump of the same size 
at the same time as the front wheels. As a further ex- 
ample of the lack of relationship between the nose wheel 
load from static distribution with that under dynamic 
conditions, move the rear wheels forward to a posi- 
tion directly under the c.g. The load on the nose wheel 
by Eq. (1) is now zero. Yet, if the airplane strikes on 
the nose wheel first, the wheel must still decelerate the 
weight, w, and the nose wheel load will be determined 
by Eq. (7): Similarly, if the airplane were to taxi in 
this condition, the nose wheel would have to lift the 
weight, w, over the bumps. Thus, design criteria based 
on static load distribution may give no indication of the 
true load on the nose wheel. 

To determine the usual magnitude of the effective 
weight, w, the radii of gyration of a number of air- 
planes were plotted against the airplane length, L, 
in Fig. 2. For single-seater, single-engined airplanes 
with the engine in the nose, it was found that the radius 
of gyration is very nearly 0.184 times the length, or 




















k = 0.184L (10) 
Substituting in Eq. (8) 
he | (0.184L)? | 
(0.184L)? + 1,? 
or 
w 0.034827 0.0348 (11) 
W  0.0348Z2 + 1,2 0.0348 + (1,/L)? 
For 


Values of w/W are plotted against 1,/Z in Fig. 3. 
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Fic. 2. Variation of pitching radius of gyration with overall 
length for conventional, single-engined military airplanes. 


the usual tricycle gear airplane, /; is between 20 and 35 
per cent of the airplane length. From Fig. 3 this means 
that between 46 and 22 per cent of the airplane weight 
acts on the nose wheel. The static percentage is 
usually only about 15 per cent, so that the error in 
designing from the load in the three-wheel attitude may 
exceed 200 per cent, assuming the same sinking speeds 
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Fic. 3. Variation of percentage airplane weight acting on the 
nose wheel with percentage overall length of moment arm from 
nose wheel to center of gravity (for conventional, single-engined 
military airplanes only). 
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in both conditions. It should be remembered that 
these figures apply only to single-engined pursuit air 
planes with the engine in the nose. Data on other air- 
planes indicate that the relationship between the two 
types of loading varies widely with the airplane type, 
purpose, loading condition, and size. The ratios be 
tween nose first and three-wheel percentages of weight 
acting on the nose wheel have been found to vary from 
4:1 to 1'/2:1. 

To relate the single-engined tricycle gear airplane 
to something more familiar, compare it with the aver 
age automobile.? Such an airplane, with a wheel base 
of 35 per cent of its length and with 14 per cent of its 
weight carried in static position by the nose wheel, 
bears the same relation to its wheel base as a conven- 
tional automobile with a normal 120-in. wheel base 
would have to a 4-ft. wheel base with the c.g. 7 in. 
forward of the rear axle. It is easily visualized with this 
overhanging arrangement how the front wheels would 
have to lift a large percentage of the car weight in order 
to pass over a bump. 

It is possible from the equations developed to deter- 
mine the distance, /;, to give the minimum structural 
weight. Roughly it may be said that the most efficient 
position is that where the nose wheel load is the same 
for the nose first landing as for the three-wheel landing. 
Assuming the same sinking speed, this could be deter- 
mined for any particular airplane by setting Eq. (1) 
equal to Eq. (9) and solving for /;. However, it has 
already been stated that the ordinary tricycle gear nose 
wheel supports about 15 per cent of the airplane 
weight. Examination of Fig. 3 discloses that for 15 
per cent of the weight to be effective on the nose wheel 
in a nose wheel first landing, /; would be nearly 45 per 
cent of the airplane length. To be consistent with the 
assumed static distribution, the wheel base would then 
be [100/(100 — 15)]45 = 53 per cent of the length of 
the airplane. This is obviously too long to be prac- 
tical. This example, it must be remembered, applies 
only to conventional, single-engined airplanes and 
again includes the doubtful assumption that the 
vertical velocity is the same for both types of landings. 
If a smaller sinking speed is used in the nose wheel 
landing or on types with a more compact mass dis- 
tribution, it may be possible to reach the optimum 
condition. 

It is obvious, however, from Fig. 3, that the nose 
wheel should be put as far forward as the geometry of 
the airplane will permit. The increase in load with 
decrease in /; is so rapid as to impose a severe penalty 
on short-nosed airplanes. Because of this penalty it is 
usually wise to make the travel of the nose wheel con- 
siderably greater than that of the main wheels. If 
this is done, the acceleration in a nose first landing 
will be reduced, decreasing the loads on the nose wheel 
and structure in proportion. The nose wheel and tire 
should be selected to have a static rating based on the 
effective weight, w. Consideration should be given 
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to the ground contact area of the tire chosen. It is 
desirable that the ground pressure under the effective 
weight should not be any greater than that of the 
main tires, otherwise, there will be considerable danger 
of the tire digging in and burying itself in soft ground 
with resultant high-bending moments on the nose strut. 
The fuselage structure must also be designed to take 
the nose first landing load. This load will usually only 
be important near the nose wheel, since the high inertia 
forces of the items in the nose will reduce the loads to 
be carried aft. 

Considering the other two items affected by nose 
wheel location—the stability against overturning and 
resistance to porpoising—both are improved by moving 
the nose wheel forward. The effect on turnover stabil- 
ity can be found under the subject of ‘‘Wheel Tread”’ 
in this paper (Part 1). The theory of porpoising has 
been given by Wylie.* The subject will also be briefly 
discussed under the section on ‘“‘Rear Wheels Location” 
in this paper (Part 1). 

One other item should be mentioned as limiting the 
nose wheel position in single-engined tractor airplanes— 
propeller clearance. The nose wheel cannot be placed 
where strut deflections or tire failure may permit it to 
enter the propeller disc. Consideration must also be 
given to providing adequate propeller ground clearance 
with any possible combination of flat tires. 

All the information assembled herein indicates that 
it is not only desirable but almost essential that the 
nose wheel be located as far forward as permitted by 
the geometry of the airplane. The only possible un- 
desirable result of this is propeller interference on single- 
engined tractor airplanes. All other characteristics 
will be improved within the practical range of nose 
wheel positions. The only exceptions to this rule 
would be airplanes with small pitching radii of gyration, 
where the nose wheel should be located so as to produce 
the same load in nose first and three-wheel landings. 


Rear Wheel Location 


The distance of the main wheels aft of the center of 
gravity is one of the most important dimensions of a 
tricycle landing gear because it affects more of the 
properties than any other single characteristic. The 
items known to be influenced by the position of the 
rear wheels are the directional stability, the longitudi- 
nal stability, the ability to land with excess air speed 
and high-sinking speeds without bouncing, the resist- 
ance to overturning, the length of take-off run, and the 
static reaction on the nose wheel. In practice, this di- 
mension is usually confined to a narrow range by two 
conflicting requirements. If the airplane is to land 
with excess air speed without bouncing, the wheels 
must be behind the airplane c.g. when they and the 
tail buffer are against the ground. Then, in any land- 
ing attitude between a tail first and a three-wheel 
landing, the vertical forces acting at the rear wheels 
will tend to decrease the angle of attack of the airplane, 
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thus reducing the lift and preventing it from rising 
into the air again. 

Frequently, it is desirable to make take-offs with the 
nose wheel in the air, thus increasing the percentage of 
air-borne weight and reducing the wheel drag. To 
do this satisfactorily the airplane must be stable longi- 
tudinally when running on the rear wheels, otherwise, 
it would either run on all three wheels or on the tail 
buffer and the rear wheels. To be stable, the condition 
is that any change in the angle of attack must produce 
a moment on the airplane tending to restore it to its 
original attitude. Moments acting on the airplane 
while taxiing result from the force of gravity, air forces, 
or wheel friction. Consider first the effect of only the 
gravity force on the airplane. This might be con- 
sidered the determination of the stability at zero air 
speed. Assume the nose wheel lifted off the ground by 
a down load on the tail of just sufficient magnitude to 
balance the moment of the airplane about the rear 
wheels. This down load is best considered as a weight 
applied to the tail as in Fig. 4 in lieu of the air load 


AIRPLANE C.G. 

RESULTANTS OF AIRPLANE 
WT. & TAIL LOAD 

Ada 










that would be obtained. by forward velocity. The re- 
sultant of the airplane and tail load then passes through 
the wheel axle and is balanced by the wheel load, P,,. 
If, however, the angle of attack is increased by a small 
angle, da, the resultant load moves aft of the wheel 
axle a distance hda and produces a moment 


dM, = Pyhda 
or 


adM,/da = Pyh (12) 


This change in moment tends to increase the angle of 
attack still further. Therefore, the gravity force is 
destabilizing or the airplane is unstable in this condi- 
tion at zero air speed. 

Next, consider the airplane running, balanced on 
its rear wheels only as shown in Fig. 5, at an air speed, 
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Conditions when taxiing on the rear wheels. 
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V, sufficient to produce a dynamic pressure, g. It is 
assumed that the wheel friction is equal to the wheel 
load, Py, times a friction coefficient, 1; that the thrust 
is equal to the sum of the wheel friction and aero- 
dynamic drag; and that the thrust and: aerodynamic 
drag produce negligible net moments about the c.g. 
The rate of change of wing-lift coefficient with airplane 
angle of attack, dC,/da, includes the ground effect; 
the rate of change of tail lift coefficient with airplane 
angle of attack, dC,7/da, includes the effect of wing 
down-wash and ground effect. 

A change in angle of attack, da, will produce a change 
in lift on the wing (wing area = S) 


dP; = Sq(dC,/da)da 


or 


dP;/da = SqdC,/da (13) 


Similarly, there will be a change in lift on the tail 
(tail area = Sr) 


dP yr = S7rq(dCir ‘da)da 
or 


dP r ‘da = SrgqdC, 7/da 
These two increases in the lift on the airplane reduce 


the load on the wheels by an amount dP, = 
—(dPs + dPr). 


GP; @P, dP, ( dC, =) 
i tt q ra tt da Ge 


The friction on the wheel is reduced also by an 
amount dF, = udPy, = —u(dPs + dP7) or 


dF, dC, , Cyr 

ie 7 os dq * set) 
The change in moment on the airplane resulting from 
the air loads only can now be found by taking moments 
about any convenient point. The point usually used 
for checking the stability in flight is the c.g. so that 
will be used here. Accordingly, the change in moment 
due to air loads is (Fig. 5) 


dM, = cdP s = (p ih cdr, i (Lr 





(15) 


+p-—c)xX 
dP, —hdP, 


Substituting —(dP; + dP) fordP, 


dM, = cdPs+(p—c)(dPs+dPr) — (Lr +p—0) X 
dP se a hdP 

This reduces to 
aM, = pdPs = LrdP yr _ hdP (16) 


This equation, it will be noticed, could have been 
derived much more easily by taking moments about 
point A, directly above the rear wheel axles and at the 
same height as the c.g. Thus, from Fig. 5, 


dM, = pdPs — LrdPy — hdP, 


which is identical with Eq. (16). 
This illustrates an important difference between the 
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stability in free flight and the stability when running 
on the rear wheels. In free flight any external moment 
causes the airplane to rotate about the c.g. The condi- 
tion for stability is therefore that the change in air 
force moments about the c.g. for a change in angle of 
attack should tend to restore the airplane to its original 
attitude. When taxiing on the rear wheels, an external 
moment tends to rotate the airplane about the point 
A, since, for small angles, the wheels prevent A from 
moving vertically while the airplane inertia restrains it 
fore and aft. It is generally accepted that moving the 
c.g. aft decreases the stability of an airplane in free 
flight. This may be stated that as the center of rotation 
moves aft the stability of the airplane decreases and 
may even become negative. Therefore, when taxiing 
on the rear wheels, the stability is similar to what it 
would be in free flight with the c.g. moved back even 
with the rear wheel axles. As the rear wheels are 
always aft of the c.g., the stability is less when taxiing 
than in flight, and it is easily possible that an airplane 
that is stable in flight may be unstable in the rear wheel 
taxiing condition. 

The total change in moment on the airplane for a 
change in angle of attack, da, can now be found by 
adding Eqs. (12) and (16), thus combining the mo- 
ments due to gravity forces with those due to air 
forces. 


dM, + dM, = dM = Pyhda + pdPs — 
L,dP,, — hdP, 








or 
dM dPs dP ry dP, . 
nies | is h — — J = | — (17) 
ae ae ae 
Substituting Eqs. (13), (14), and (15) in Eq. (17) 
dM dC, dC, 
—- = P,)} S — LS — 
da bal * “Sith acl tied 
wi(s8 4 5,8e2)] as 
da da 


For the airplane to be stable, dM//da must be nega- 
tive—that is, an increase in angle of attack must 
produce a negative moment that will tend to return the 
airplane to its original attitude. An examination of 
Eq. (18) shows that the first term, P,,h (the load on 
the wheels times the height, 1), is always positive, so 
that for stability the term in the brackets must be 
negative. If this term is negative, then the stability 
improves with g or in proportion to the velocity 
squared. At low air speeds, the airplane is always 
unstable for the condition of running only on the rear 
wheels. This is not important, however, as at low air 
speeds there is not enough tail load available to lift the 
nose wheel off the ground. 

In practice, it is usually satisfactory if the airplane 
becomes stable before it reaches the take-off speed. 
As the airplane approaches the take-off speed, the load 
on the wheels, P,,, becomes smaller and smaller, so 
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that at the instant of take-off the term P,,h in Eq. (18) 
becomes zero. Therefore, to be stable at take-off the 
condition is that the term in the brackets must be 
equal to, or less than, zero. 

















SdC, dCir ( dC, Hae) “ 
_ — j = 0 
P da Lets da 7 ~~ ¥ & da 
or 
d dC. dC dC, 
bet he Hn (si + Sr ss) (19) 
da da da da 


If Eq. (19) is not satisfied, the airplane will become 
progressively less stable as the speed increases. This 
would necessitate the pilot holding the nose wheel on 
the ground until take-off speed was reached, at which 
any slight disturbance would produce a sudden in- 
crease in the angle of attack, causing the airplane to leap 
into the air. If sufficiently unstable, it might be im- 
possible for the pilot to hold the nose on the ground 
until take-off speed was reached. Should the pilot at- 
tempt a nose high take-off with an unstable airplane, 
the angle of attack would increase until the wing 
stalled. Beyond the stall, the slope of the wing-lift 
curve becomes negative so that the terms involving 
dC,/da become negative, the equation is satisfied, and 
the airplane becomes stable again. In practice, the tail 
buffer usually strikes the ground and drags before the 
stall is passed sufficiently to reverse the lift curve slope. 
The airplane may then run along on the rear wheels 
and tail buffer or, if sufficiently unstable, the angle of 
attack may increase still further, leaving only the tail 
buffer on the ground. In either case, the take-off run 
would be increased by the high drag of the stalled wing 
plus that of the tail buffer, and the airplane would be 
forced off the ground in a stalled attitude. For the 
case of take-off on the tail buffer alone (which can occur 
only on airplanes where the buffer is quite well for- 
ward), the airplane would become laterally unstable 
as soon as the rear wheels left the ground and would 
be very likely to go up on one wing tip. Thus, it can be 
seen that while a tricycle gear airplane so constructed 
as to be unstable with the nose wheel in the air could be 
flown it would be a highly treacherous and unsatis- 
factory arrangement. 

Eq. (19) is not in a form that is easily used, since it 
involves the rates of change of lift coefficients with air- 
plane angle of attack. As these are usually unknown, 
the formula will be transformed into a form using the 
lift curve slopes of the airfoils for infinite aspect 
ratio. 


Let 

a = slope of wing-lift curve = dC,/da for infinite 
aspect ratio 

a7 = slope of tail lift curve = dC,7/da for infinite 
aspect ratio 

R = wing aspect ratio 


Rr = tail aspect ratio 


The slope of the lift curve of a wing of aspect ratio, 
R, in terms of the slope for infinite aspect ratio 
is! 

dC,/da = a/\1 + [a(1 — o)/xR)}} (20) 


where o is a factor depending on the ratio of wing 
height above the ground to wing span and is available 
in standard texts on aerodynamics. 

The tail operates in the down-wash of the wing. 
The change in the tail angle of attack for a given change 
in airplane angle of attack is smaller by the change in 
the down-wash angle, d8. Or 


darait = deg irptane ‘03 dg (21) 


When the ground effect is considered, d8/da at the 
tail varies with its height above the ground. If the tail 
is approximately even with, or above, the wings! 


dB/da = (2a/rR)(1 — a) (22) 
Accordingly, using Eqs. (20), (21), and (22), the 


rate of change of tail-lift coefficient for airplane angle 
of attack becomes: 


dCir ae FE ar |[: Sy 2a (1 ) 
da 1+ (ar/eRz)(1 — or) rR % | 


where o7 is the ground effect factor for the tail. 
Eq. (19) may now be rewritten, using Eqs. (20) and 
22), as 








Srar [1 + (a mR)(1 _ o)] 


Ss 11+ @ec/eRD — esl * 
97 — 
Pecan) ere 


n(Lr — uh) 


where 7 is a tail efficiency factor to account for inter- 
ference effects. 

The stability of the airplane may be determined by 
Eq. (24), since all of the quantities are known except 7, 
which may either be determined from wind-tunnel 
tests or estimated from experience. However, for con- 
ventional airplanes this equation may be simplified 
considerably by making a few slightly conservative as- 
sumptions. The wing height above the ground is 
usually between one-tenth and two-tenths of the wing 
span. Within this range' o has values between 0.655 
and 0.485. The effect of o on the final result is small, 
and it will usually be slightly conservative to assume ¢ 
equal to 0.5. The height of the tail surface above the 
ground varies with the angle of attack of the airplane, 
but in level position the tail is usually about one-half 
its span above the ground. For this attitude o7 may 
be assumed zero with only a small conservative 
error. 

As the attitude of the airplane changes, moving the 
tail closer to the ground, o 7 will increase, thus improv- 
ing the stability. The maximum possible improve- 
ment, if the tail were just above the ground line, 
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would be 200/R,r per cent, which, for usual tail aspect 
ratios of about four, would be not more than 50 per 
cent. In most cases, ¢7 = 0 will be a satisfactory ap- 
proximation for the most severe condition with the 
tail at maximum height. The slope of the lift curves 
of both the wing and tail for infinite aspect ratio may 
be assumed to be equal—a condition that is very nearly 
true for modern airfoils. In addition, they may be set 
equal to the theoretic value of 27. This value is 
slightly higher than usually realized, but it can be seen 
from Eq. (24) that the actual numerical value is of 
secondary importance so long as the correct ratio is 
kept between the values for the wing and tail. Two 
other items—z , the coefficient of friction, and 7, the 
tail efficiency—can be replaced by constants. The co- 
efficient of rolling friction? on airports in fair condition 
is about 0.03. 

Tail efficiencies on modern airplanes are usually 
70 to 80 per cent; so for preliminary design around 
work a value of 0.7 should be satisfactory. Substi- 
tution of these terms in Eq. (24) gives: 


[1 + (2r/xrR)(1 — 0.5)] 
[1 + (24/rR_)(1 — 0)] 
E = sot = 0-8 > p + 0.03h 


us 





S 
0.7 4 (Ly — 0.03h) 


This reduces to 


Rr(R — 1) 


—_——— — 0.034 = 
R(Rr + 2) =? 


S 
0.7=5 Lr — 0.03h) (25) 


The term (L; — 0.03h) is the effective tail arm. 
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Graph for determination of static longitudinal stability 
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when taxiing on the rear wheels. 


The second part of this term is small compared to L, 
and can be ignored. 

Also, if a tail effectiveness factor, e, is substituted for 
R7(R — 1)/R(Rr + 2), Eq. (25) becomes 


(0.7eS 7 S) - (0.03h Lr) Pa pb Lr 


In Fig. 6, values of € are plotted against wing aspect 
ratio, R, for constant values of tail aspect ratio, Rr. 
Fig. 7 has values of p/L, plotted against «S;/S for 
constant values of h/L,y. The entire solution for the 
maximum value of » may thus be done graphically 
if only the wing aspect ratio and area, tail aspect 
ratio and area, tail length, and c.g. height are known. 
The answer as finally determined should be satisfac- 
tory for design work in conventional airplanes. If the 
airplane lies outside the range of assumptions, the stabil- 
ity should be determined directly from Eq. (24). It 
should be remembered in applying this equation that 
the distances p and Ly are always measured parallel 
to the ground line, while / is measured perpendicular 
to it, regardless of the attitude of the airplane. In 
practice, 17 and / are little affected by angle of attack, 
but p usually decreases rapidly as the angle of attack 
increases. Therefore, as the angle of attack is in- 
creased, the stability improves. 

By making the same substitutions and approxima- 
tions in Eq. (18) as were made on Eq. (19) and by re- 
placing the wheel load, Py, by the airplane weight ‘ess 
the wing lift, or (W — gSC,), the dynamic pressure, 4, 
and consequently the speed at which the airplane be- 
comes stable may be calculated for any particular 
angle of attack. 


(26) 
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In using this equation, the values of » and C, must 
be for the same angle of attack and the ground effect 
must be considered in determining C,. 

The minimum speed at which the elevators and 
wing can lift the nose wheel into the air may be deter- 
mined if the lift coefficient of the tail is known for the 
elevator full-up position and airplane level angle of at- 
tack. At speeds lower than the minimum at which 
enough down load on the tail can be developed to lift 
the nose wheel into the air, the stability when running 
on the rear wheels only has no real meaning unless it is 
desired to consider the extreme case when the nose 
wheel strikes an obstruction large enough to throw it 
into the air. It would be desirable to have the air- 
plane stable at all possible angles of attack at this 
minimum speed, but this condition usually cannot be 
realized. The next best thing would be to have the air- 
plane stable at this minimum speed and the angle of 
attack corresponding to the attitude with the tail 
buffer just above the ground. Attainment of this con- 
dition may also be impossible, since any changes made 
to improve the stability, such as reducing p or increas- 
ing the tail arm or area, will reduce the speed at which 
the nose wheel can be lifted and thereby require a 
further stability increase. 

Little information is available on the degree of 
stability actually used on tricycle gear airplanes. 
Approximate checks on several airplanes indicate 
minimum stable speeds of the order of two-thirds to 
nine-tenths of the take-off speed. This apparently 
is not usually the result of design for stability but rather 
is the result of a desire to make it as easy as possible for 
the pilot to lift the nose into the air and to reduce the 
static nose wheel load. 

At present, it is felt that Eq. (26) provides the best 
criterion of stability. If this equation is satisfied for 
the level attitude of the airplane, the stable speed will 
be reached before the take-off speed and appreciably 
so at the higher angles of attack. No dangerous condi- 
tions should arise from instability below the take-off 
speed, the worst possibility being a banging and 
dragging of the tail buffer on the ground and a conse- 
quent increase in length of take-off run. 

The analytic treatment of the stability problem 
given herein is a modification of the method of Wylie? 
with the addition that the effect of c.g. height is con- 
sidered. As given here, the criterion developed is a 
measure of the static stability only. It may be ap- 
plied to the dynamic case if desired, but it is felt that 
there is little to be gained by such a procedure. 


) 
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The two conditions—(1) that the c.g. must be for- 
ward of the rear wheels at the maximum ground 
angle and (2) that the airplane must be stable on the 
rear wheels at take-off speed—usually limit the satis- 
factory rear wheel location to a narrow range. The 
effect of changing the wheel location within this region 
on other characteristics of the airplane will next be 
considered. . 

The directional stability is proportional to the dis- 
tance of the rear wheels aft of the c.g. The relationship 
between this distance and the directional stability is 
developed in the section of this paper on ‘Wheel 
Tread” (Part 1), where it is found that the minimum 
acceptable distance from the wheels to the c.g. usually 
gives more than adequate directional stability. 

The resistance to porpoising or longitudinal oscilla- 
tion when running on all three wheels increases with 
rear wheel distance from the c.g. This type of oscilla- 
tion usually only occurs with extremely short wheel 
bases and high c.g. positions. Generally, it appears 
only at low speeds, since the damping effect of the 
tail prevents oscillation at high speeds. The only 
detrimental effect, therefore, is possible discomfort of 
the crew. This condition is not considered important 
enough to warrant further discussion. However, if 
desired, the approximate period of oscillation may be 
calculated by the method outlined by Wylie.’ 

Resistance to overturning is improved by keeping the 
rear wheels close to the c.g. The effect is usually small, 
the most important factors in preventing overturning 
being the tread and distance from c.g. to nose wheel. 
This problem is considered in the section of this paper 
on ‘Wheel Tread” (Part 1). 

From the previous considerations, the best position 
for the rear wheels is believed to be just aft of the rear- 
most c.g. with the tail buffer on the ground. Placing 
the wheels further aft of the c.g. can only improve the 
porpoising characteristics, which are unimportant, and 
the directional stability, which is usually adequate and 
may thus become excessive. At the same time, the 
longitudinal stability during nose-high take-offs and 
the resistance to overturning will be decreased. Thus 
there is one partially desirable effect against two 
and possibly three undesirable effects. 
cept in the case when it is necessary to improve the 
porpoising characteristics at the expense of the other 
items, the rear wheels should be located as far forward 
as permitted by the rearmost c.g. 


Therefore, ex- 


WHEEL TREAD 

General 

The wheel tread is a dimension that, in single-engined 
types, is not usually limited or closely controlled by the 
airplane geometry. In addition, the characteristics of a 
tricycle gear airplane affected by the tread are also a 
function of wheel base and fore-and-aft location of the 
c.g. The determination of the wheel positions fore and 
aft has been discussed in the section of this paper on 
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‘‘Wheel Base” (Part 1), where it is found that several 
of the airplane characteristics determine the longi- 
tudinal wheel locations within narrow limits and that 
these characteristics are independent of the tread. 
Therefore, in the design of a tricycle landing gear, the 
fore-and-aft wheel positions should first be determined 
by the aforementioned criteria. The wheel tread can 
then be adjusted to make satisfactory the properties of 
the landing gear affected by both wheel base and tread. 
In multiengined airplanes where the tread is deter- 
mined by the nacelle spacing, it may be necessary to 
modify this method. The tsual nacelle spacing, how- 
ever, is sufficient to give an entirely adequate tread. 

The qualities influenced by the wheel tread are the 
resistance to overturning and the ease of steering by 
differential braking. In addition, there may be second- 
ary effects on the directional stability. 


Stability Against Overturning 


An airplane may be said to be stable against over- 
turning if the wheels will slide before they can develop 
a frictional force sufficient to upset the airplane. Con- 
sider the triangle formed by the three wheels as shown 
in Fig. 8. Let the resultant inertia and gravity force 
acting at the c.g. be represented by a vector. If this 
resultant vector intersects the ground inside the tri- 
angle, then the airplane will not overturn; if it inter- 
sects the ground outside the triangle, then it will over- 
turn. 

If 


b = tread 
P., = total wheel load (all three wheels) 
hy height of c.g. above ground line 








Fic. 8. 


The largest possible magnitude of the horizontal com- 
ponent of the force acting on the airplane cannot exceed 
uP,, the maximum frictional force of the wheels. Ob- 
viously, this force is most likely to upset the airplane 
when it acts in a direction perpendicular to one side of 
the triangle. The overturning moment is then: 


Mo = uP | (28) 
The stabilizing or resisting moment is 
M, = P»X (29) 
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If the airplane is not to overturn, M@, = Mb, or 


Pak SS ew oh OC 2A. Zz uk (30) 

But 
e ' b/2 
X =i, sin y and sin y = VP + 6/2) 

So that 

X = hb/2V 2 + (6/2)? (31) 
Substituting Eq. (31) in Eq. (30) and solving for } 

b > Quhl/Vh? — phy? (32) 
And 

b/l = Qyh,/lV (h/1)? — u2(hy/l)2 (33) 


Minimum values of tread to prevent overturning 
can be calculated from Eq. (32) if the coefficient of 
friction, wu, is known. Eq. (33) is the nondimensional 
form of Eq. (32), expressing the minimum satisfactory 
ratio of tread to wheel base in terms of the ratios of 
c.g. height to wheel base and nose wheel arm (/;) to 
wheel base. 

The coefficient of friction, u, varies with the type of 
surface and its condition. A value of 0.55 is customa- 
rily used in brake design to produce a brake capable of 
skidding the tire. A coefficient of friction of 1.00 is 
usually assumed in ground-loop conditions for design 


of landing gears. These conditions, however, are: 


usually specified so as to give certain loads on the 
structure and not necessarily to simulate actual con- 
ditions. In view of these wide discrepancies between 
the values for the coefficient of friction, no final criteria 
can be given. A coefficient of 0.85, however, is be- 
lieved to be adequate for use in this formula. This 
has a slightly better than 1.5 safety factor over 0.55, 
the coefficient used for braked friction. 

Objections have been raised to this overturning 
criterion on the ground that so long as the nose wheel is 
free swiveling the airplane will turn sufficiently in the 
direction it was yawed to prevent overturning. To 
turn, however, the rotational inertia of the airplane 
must be overcome. Accordingly, computations were 
made of the lateral acceleration of airplanes landing 
with various angles of yaw by assuming the side 
load resisted only at the rear wheels. It was found 
that angles of yaw above approximately 10 deg. would 
be sufficient to cause the rear wheels to skid. The re- 
sisting frictional force would, of course, be proportional 
to the vertical load on the rear wheels rather than on all 
three wheels. Since the rear wheels usually support 
about 85 per cent of the airplane weight, the force 
necessary to slide the wheels is about 85 per cent of 
what it would be to slide all three wheels. However, 
as mentioned previously, the maximum value of the 
coefficient of friction is about 1.00. Using 0.85 as was 
done in computing the curves, will, therefore, approxi- 
mately compensate for the swiveling of the nose wheel. 
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Minimum ratio of tread to wheel base which will prevent 
overturning for a coefficient of friction of 0.85. 
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Values of 6/1 have been plotted against /,// for 
constant values of /,// in Fig. 9, using » = 0.85. These 
curves may be used for other values of u by multiplying 
the values of h,// by u/0.85 and getting an equivalent 
h,/l. Read the value of b/] from the curves using the 
equivalent value of h,/1. 

The curves are labeled in terms of /;//, where J, = 
(1 — h), instead of 1,/] as indicated by the formula, 
since /,/] corresponds to the percentage of static load 
on the nose wheel and is usually more convenient than 


1,/l. 


Directional Stability and Ease of Steering 


The directional stability is determined by the dis- 
tance of the rear wheels aft of the c.g. and the height 
of the c.g. above the ground line. Neither of these di- 
mensions is related to the tread. Ease of steering may 
be said, however, to be inversely proportional to the 
directional stability and, if the steering is done by 
differential braking, directly proportional to the tread. 
It is difficult to give a criterion for directional stability 
alone because, since all tricycle landing gears are di- 
rectionally stable so long as the rear wheels are behind 
the c.g., the problem becomes one of the degree of 
directional stability desired. It has been pointed out 
in the section of this paper on wheel base (Part 1) 
that the value of /, is usually controlled by requirements 
distinct from those of directional stability. Similarly, 
h, is usually determined by the airplane geometry and 
cannot be altered to adjust the directional character- 
istics. The directional stability is, therefore, usually 
set by other requirements of the airplane, and this de- 
gree of stability must be accepted as a characteristic 
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of the airplane. This is not so serious as it might at 
first seem, since airplanes can be operated satisfactorily 
over a wide range of directional characteristics. The 
conventional landing gear is unstable until the air 
speed becomes sufficient for the stabilizing character- 
istics of the tail to 6vercome the instability of the 
landing gear. Yet little difficulty is experienced by 
pilots in overcoming ground-looping tendencies except 
in airplanes with extremely narrow treads, where the 
wheel brakes cannot supply an adequate stabilizing 
moment, or under other adverse conditions. As the 
tricycle landing gear is always stable, ground looping 
cannot occur, and it is never necessary for the pilot to 
supply a stabilizing moment. Instead, the problem is 
reversed; in order to make a turn the pilot must pro- 
vide a moment to partially overcome the naturally 
stable characteristics of the airplane. Accordingly, the 
sharpness of turn depends on the magnitude of the 
turning moment that can be applied which, assuming 
the turning is done by differential braking, depends 
on the wheel tread and tire coefficient of friction. 

Consider a tricycle gear airplane of mass, W/g, 
following a path of radius of curvature, R, at a velocity, 
V, as shown in Fig. 10. There will be a centrifugal 
force acting outward at the c.g. of 


F = (W/g)V?/R (34) 
—_—" Re dee 
a l 5 





Fic. 10. 


The nose wheel is free to swivel and, therefore, can- 
not resist any of the. side force, F, which accordingly 
must be taken entirely by the rear wheels. This pro- 
duces a moment tending to return the airplane to a 
straight path, of 


M, = hF = 1.WV?/gR (35) 


The force, F, also acts at a height, 4, above the 
ground line, thereby producing a moment tending to 
increase the load on the outside wheel and reduce it on 
the inside. The change in load is (using Eq. (34)) 


AP = Fh,/b = h,WV?/bgR (36) 


This change in load will increase the wheel friction 
on the outside of the turn and, accordingly, produce 
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another moment tending to restore the airplane to a 
straight course. The reduction of load on the inside of 
the turn will also have the same effect. The airplane, 
following the path of radius, R, therefore has a total 
moment tending to restore it to a straight course equal 
to Eq. (35) plus Eq. (36), which reduces to 


M, = (WV?2/gR)(l, + uli) 


Ii steering is done by differential braking, the 
moment produced by braking the inside wheel is 


(37) 


Mp = (6/2)(uz — wo)Pr (38) 


where yz is the coefficient of friction of the braked wheel 
as distinguished from yo, the coefficient of friction of the 
unbraked wheel, and P; is the load on the inside wheel. 
But P; is equal to the static load on the wheel less 


AP. 


P, = [(W/2)h/l] — [mWV?2/bgR] (39) 


So that substituting Eq. (39) in Eq. (38) 


(b/2)(ue — mo)W{h/2l — (ii V?/bgR)| (40) 


Ms = 
To make a steady turn, M/, must equal Ws, or using 
Eq. (37) and Eq. (40) 


(WV2/gR)(l2 + wo) = (6/2)(uz — wo) X 


W[h/2l — (i, V2/bgR)] 


Solving for R 





lp + '/ohi(up + | (41) 


Mp — Bo 


R = [4V?l gon | 


The value of R may now be determined by substitut- 
ing numerical values in Eq. (40). The sharpest turn 
will, of course, be made by sliding the inside wheel. 
For this case ug will be assumed to be equal to 0.55 
and po to 0.03. The equation then becomes 


v2] 
R = 0.514 (lp + 0.29h1) (42) 


bh, 


where V is in m.p.h. and the rest of the dimensions are 
in ft. 

As an example, consider a typical pursuit airplane, 
with the following dimensions, taxiing at 30 m.p.h.— 
b = 8.33 ft., 2 = 10 ft., 4, = 8.75 ft., , = 1.25 ft. and 
h, = 5.0 ft. Then R = 172 ft. 

Thus, in this case, the airplane is so stable direc- 
tionally that at 30 m.p.h. with one main wheel locked 
it takes a 350-ft. circle to turn around. The diameter 
of this circle is proportional to the square of the speed 
so that the size for any other speed can be found by a 
simple ratio. At a stalling speed of approximately 
80 m.p.h., the diameter would be nearly 2,500 ft. Thus, 
no serious consequences would result from locking one 
wheel, even when landing, since the turning radius is 
so large. It is impossible to tip a tricycle gear airplane 
up on one wing by differential braking because, as the 
radius of turn becomes smaller, the load on the inside 
tire which is providing the turning moment decreases, 
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thus reducing the frictional force and the turning 
moment. For the airplane to go up on one wing, the 
inside wheel would have to lift clear of the ground, 
reducing its turning moment to zero. The frictional 
drag of the wheel on the other side would then pro- 
duce a turning moment in the opposite direction which 
would result in a centrifugal force, forcing the braked 
wheel back onto the ground. It is possible, however, 
by increasing the tread, 0, sufficiently to make R so 
small that the centrifugal force is greater than the 
resistance of the tires to sliding sidewise. That is, if 
as is the coefficient of friction for sliding sidewise, 


(W/g)V2/R = u.Wh/l 


Or, to prevent sliding 


R= V?l/usgh (43) 


If u, is taken to equal 0.55, the airplane of the example 
would begin to skid at 30 m.p.h. for a radius of turn of 
125 ft. Thus, no matter what the tread was increased 
to, the airplane could not turn around in a circle less 
than 250 ft. in diameter at 30 m.p.h. Since the turn- 
ing radius is inversely proportional to the tread, the 
corresponding value of tread required would be 6 = 
172/125 XK 8.33 = 13.75 ft. The same result could 
have been obtained by setting Eq. (43) equal to Eq. 
(41) and solving for 0. 

In applying these equations, it should be remem- 
bered that, in their development, R was considered 
large in relation to the tread. It can be seen from 
Fig. 10 that for small values of R the frictional drag of 
the nose wheel would have a stabilizing effect and the 
centrifugal force, F, would not be perpendicular to the 
planes of the rear wheels. For radii greater than five 
times the tread, these discrepancies have a negligible 


effect. 
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